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PREFACE 


As  a part  of  the  investigation  of  a possible  Sea  Level  Canal  con- 
necting the  Atlantic  and  Pacific  Oceans,  the  U.  S.  Army  Engineer  District, 
Jacksonville,  requested  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  to  investigate  roughness  that  would  result  in  the  canal. 

This  report  was  prepared  by  Dr.  G.  H.  Keulegan  with  the  assistance 
of  CPT  M.  J.  Mathews,  formerly  with  the  Mathematical  Hydraulics  Division 
of  the  Hydraulics  Laboratory,  WES. 

The  work  on  this  study  was  commenced  in  1969  and  was  followed 
somewhat  intermittently  since  then.  The  study  was  under  the  general 
supervision  of  Mr.  E.  P.  Fortson  (retired),  former  Chief  of  Hydrau- 
lics Laboratory,  and  Mr.  H.  B.  Simmons,  present  Chief;  Mr.  F.  A. 

Herrmann,  Jr.,  former  Chief  of  Estuaries  Division,  and  Mr.  R.  A.  Sager, 
present  Chief. 

Directors  of  WES  during  the  course  of  this  investigation  and  the 
preparation  and  publication  of  this  report  were  COL  J.  R.  Oswalt,  Jr., 

CE,  COL  L.  A.  Brown,  CE,  BG  E.  D.  Peixotto,  CE,  COL  G.  H.  Hilt,  CE, 
and  COL  John  L.  Cannon,  CE.  Technical  Directors  were  Messrs.  J.  B. 

Tiffany  and  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con' 
verted  to  metric  (SI)  units  as  follows: 


Multipl; 


To  Obtain 


0.02831685 
0. 017^5329 
0.301+8 
0.301+8 
0.301+8 
25.1+ 

1.60931+1+ 

0. 1+535924 


cubic  feet 
degrees  (angle) 
feet 

feet  per  second 

feet  per  second  per  second 

inches 

miles  (U.  S.  statute) 

pounds  (mass)  per  second 
per  foot 

square  feet 

square  feet  per  second 


cubic  metres 


0.09290301+  square  metres 

0.09290301+  square  metres  per  second 


AN  ESTIMATE  OF  CHANNEL  ROUGHNESS  OF  INTEROCEANIC  CANALS 


PART  I:  RESISTANCE  OF  NONALLUVIAL  CHANNELS 


Introduction 


1.  The  present  study  deals  with  the  expected  values  of  the  rough- 
ness of  the  interoceanic  canals  created  by  either  nuclear  or  conventional 
construction  methods.  Nuclear-cratered  channels  pose  two  problems.  The 
first  refers  to  the  roughness  under  the  hypothetical  assumption  that  the 
crater  channels  are  of  a uniform  cross  section  throughout.  Knowing  the 
distribution  of  rock  sizes  in  the  ejecta  and  in  the  fallback,  the  rough- 
ness is  to  be  inferred  in  terms  of  Manning's  n . Eventually,  a rela- 
tion of  rock  size  and  Manning's  n can  be  established  by  relying 

upon  the  concept  of  equivalent  sand  roughness.  Next,  with  this  informa- 
tion as  a basis,  the  effective  Manning's  n of  crater  channels  is  de- 
rived subject  to  the  known  distribution  of  rock  sizes  and  uniform  cross 
section.  The  second  problem  deals  with  the  increase  of  channel  resis- 
tance resulting  from  the  expansions  and  contractions  of  a nuclear-blasted 
channel.  Part  of  the  answer  to  this  question  comes  from  the  experiments 
conducted  at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES). 
Difficulties,  however,  arise  in  the  attempt  to  apply  laboratory  results 
to  prototype  situations.  The  ordinary  model  rules  for  the  transfer  of 
roughness  do  not  apply  since  the  alternating  expansions  and  contractions 
bring  about  losses  that  are  not  related  to  wall  skin  friction  but  are 
caused  by  internal  turbulent  friction.  The  proper  transfer  relation  is 
established  by  developing  the  theory  of  losses  encountered  in  the  mixing 
of  a current  moving  over  adjacent  stagnant  or  static  waters. 

2.  In  canals  created  by  conventional  methods  the  bottom  is  com- 
posed of  Atlantic  muck.  Eden,  formerly  of  the  U.  S.  Army  Engineer 
District,  Jacksonville  (SAJ),  has  informed  us  that  the  muck  has  the 
following  grain  size  distribution:  100  percent  less  than  0.1  mm, 

65  percent  less  than  0.01  mm,  and  35  percent  less  than  0.002  mm.  This 
extremely  fine  composition  poses  the  difficult  problem  of  deciding 


exactly  what  the  equivalent  roughness  would  he.  In  ordinary  alluvial 
streams,  with  bed  materials  of  a mean  grain  size  on  the  order  of  a 
millimeter  or  a tenth  of  a millimeter,  the  resistance  is  related  to  the 
presence  of  ripples  or  dunes.  For  flows  with  a Froude  number  less  than 
0.5,  and  if  the  velocity,  grain-size  distribution,  channel  dimensions, 
and  channel  slope  are  known,  the  corresponding  resistance  can  be  roughly 
established.  The  question  of  bed  feature  formation  with  very  fine  sand 
has  rarely  been  investigated.  If  it  is  found  that  the  effect  of  ripples, 
if  present,  is  of  no  consequence,  the  channel  boundaries  may  be  regarded 
as  hydrodynamic ally  smooth  and  having  small  friction. 

3.  The  importance  of  the  problem  dictates  that  an  effort  be  made 
to  determine  the  condition  or  the  likelihood  of  the  formation  of  bed 
feature  with  fine  particles  and  the  resultant  resistance.  Two  ap- 
proaches are  open  for  the  examination  of  friction.  First,  various 
theories  on  the  resistance  of  alluvial  streams  with  bed  materials  of 
moderate  sizes  are  considered.  Selecting  from  these  the  one  most  con- 
venient for  the  purpose,  relations  are  projected  to  canals  with  fine 
bed  material.  Secondly,  an  attempt  is  made  to  determine  the  resistance 
of  individual  ripples  or  dunes  assuming  that  the  relevant  current  veloc- 
ity is  the  velocity  at  the  apices  of  the  ripples  or  dunes.  The  resis- 
tance coefficients  are  established  for  each  form  by  referring  to  labora- 
tory results  which  have  noted  the  dimensions  of  these  forms  in  addition 
to  the  stresses  at  the  walls  and  the  bottom.  If  there  is  additional 
information  on  the  sizes  of  these  bed  features  in  fine  sediment  channels, 
the  resistance  coefficient,  being  a dimensionless  quantity,  could  be 
applied  directly. 

Estimates  of  Roughness  of  Crater  Surfaces 

An  expose  of  open  channel  formulae 

1+.  An  earlier  study  on  the  laws  of  turbulent  flow  in  open  chan- 
nels  led  to  the  relation 


TT  = 6.25  + 5.75  log  (1.1) 


where 


(1.27 


U = channel  mean  flow  velocity* 

U*  = shear  velocity 

g = gravitational  acceleration  constant 
R = hydraulic  radius 
S = the  energy  gradient 

In  the  event  that  the  depth  is  uniform  and  the  variation  of  velocity  head 

with  distance  is  negligible,  S merely  represents  the  slope  of  the  water 

surface;  k = equivalent  sand  roughness  height.  Equation  1.1  represents 
s 

the  logarithmic  law  of  mean  velocity.  This  is  also  nearly  equivalent  to 


provided  that 


(1.3) 


(l.M 


which  is  a form  of  Manning's  law.  The  ordinary  form  of  this  law  when 
lengths  are  measured  in  feet**  is 


u = 1 • r1//6/rs 

n 


(1.5) 


which  may  be  written  also  as 


U = 


0.260 


/gRS  R 


1/6 


where  R is  expressed  in  feet.  Accordingly 


U _ 0.260  „l/6 
7?  — n 

U*  n 


(1.6) 


(1.7) 


* For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  A). 

**  A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Si)  units  is  presented  on  page  3. 


Hence  comparing  Equation  1.7  with  Equation  1.3  gives 


n = 0.0324k 


(1.8) 


a relation  between  Manning's  n and  the  equivalent  sand  roughness  kc 
the  latter  expressed  in  feet. 


Equivalent  sand  roughness 
and  asperity  concentration 

5.  The  very  fruitful  idea  of  equivalent  sand  roughness  is  due  to 
2 

Schlichting.  The  practical  implication  of  the  idea  is  as  follows.  One 
may  consider  two  channels  or  two  conduits  having  the  same  form  and  hy- 
draulic radius.  The  walls  of  one  channel  are  covered  with  asperities 
of  height  k and  the  other  with  closely  packed  sand  of  height  kg  . 

If  the  flows  in  the  two  channels  under  the  same  pressure  difference  are 
equal,  then  ks  is  the  equivalent  sand  roughness  of  the  asperities  k 
of  the  first  channel.  Schlichting  conceived  an  ingenious  apparatus  that 
eliminated  the  use  of  towing  tanks  to  investigate  the  resistance  of  bat- 
tleship armor  design.  The  apparatus  consisted  of  a wide  rectangular 
channel,  three  sides  of  which  were  smooth  surfaces  while  the  fourth  side 
consisted  of  a rough  plate.  The  rough  plate  could  be  covered  with  any 
type  of  asperities  and  concentration  values.  The  asperities  concentra- 
tion is  defined  by  Schlichting  as  follows.  If  is  the  total  value 

of  the  projected  area  of  all  the  asperities  found  on  an  area  A , upon  a 


a plane  parallel  to  the  direction  of  flow,  then  A^/A  is  to  be  taken 
as  the  concentration  of  the  asperities.  If  the  projection  of  a single 
asperity  is  a^  and  there  are  N number  of  them  on  an  area  A , then 

A = a N . For  a layer  of  spheres  of  diameter  2a,  spheres  squarely  com- 

r r 2 2 
pacted,  taking  N = 1 , A = 4a  , A^  = tra  , and  A^/A  = 0.783  • 


Schlichting  recognized  that  the  relation  between  k and  k 
the  form 


the  function  relation  being  different  for  different  forms  of  the  asperi 
ties.  For  example.  Figure  1 shows  the  manner  of  dependence  for  spheres 
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Figure  1.  Effect  of  asperities  concentration  on  resistance 

as  found  by  Schlichting.  Also  shown  in  this  figure  are  the  resistance 

characteristics  of  cubes.  Part  of  the  data  of  the  cubes  was  taken  from 
3 4 

Iowa  tests  ’ and  the  remaining  part  from  WES  tests.  The  investigators 

of  the  latter  tests.  Miller  and  Petersen,  do  not  give  the  quantities 

A /A  and  k /k  in  their  report,  but  they  do  give  the  sizes  of  the 
r s 5 

cubes  and  the  spacing  of  the  cubes.  Among  other  things  they  give  the 
depth  of  the  water  y^  , hydraulic  radius  r , and  Manning's  n . The 
values  of  depend  upon  the  values  of  k and  the  spacing.  The 

expression 

5.75  log  z-.s&s.  b1/6-6.25 


which  results  upon  elimination  of  U/U*  from  Equations  1.1  and  1.7 

enables  us  to  determine  k from  the  known  values  of  R and  n . 

s 

Forming  ks/k  and  y^/k  , the  limiting  values  of  ks/k  for  greater 
water  depths  are  determined  from  a plot  as  in  Figure  2.  Referring  again 
to  Figure  1 it  is  to  be  inferred  that  for  small  values  of  A /A  the 


T 


Figure  2.  Dependence  of  effective  roughness  on  water  depth  in  channel 

ratio  ^s/k  increases  with  Ar/A  . The  ratio  kg/k  attains  its  great- 
est value  at  about  Ar/A  = O.U  . Beyond  that  there  is  a decrease  with 
increasing  A^/A  . 

Equivalent  sand  roughness 
of  rocks  and  pebbles 

6.  For  gravelly  and  rocky  channels,  data  similar  to  those  shown 
in  Figure  1 are  needed.  This  information  is  shown  in  Figure  3 for 

k 6 2 

sand,  limestone  rock,  and  Hanover  sand.  Schlichting' s experi- 
ments were  performed  on  Hanover  sand  similar  in  compactness  to  that  of 
Nikuradze. 

7.  The  value  of  ks/k  found  by  Schlichting  is  1.6  which  suggests 
that  all  sands  are  not  alike  with  respect  to  their  resistance  character- 
istics. The  full  curve  shows  mainly  the  effect  of  concentration  on  the 

value  of  k /k  . Thus  for  rocks  and  pebbles  with  concentrations  A /A 
s r 


- :r  * '~-~  r * • 
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► # 
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greater  than  0.5»  kg/k  is  nearly  constant.  When  A^/A  is  greater 
than  0.5  we  shall  regard  the  rocks  as  closely  packed  and  assume  that 
kg/k  is  1.75,  i.e.,  the  equivalent  grain  size  is  1.75  times  the  grain 
size  of  rocks. 

The  Manning's  n of 
rocks  of  a specified  size 

8.  Since  kg/k  equals  1.75  for  closely  packed  rocks.  Equa- 
tion 1.8  becomes 

n = 0.032Ml.75k)1/6 
or 

n = 0.0355kl/6  (1.: 

where  k is  the  grain  size  of  rocks  expressed  in  feet.  This  relation 
is  shown  in  Figure  b.  The  upper  curve  is  for  the  range  0.1  < k < 10  ; 
the  lower  curve  is  for  the  range  0.001  < k < 0.1  . Thus  a 

10 


•L, 


• . . i 


Figure  Relation  between  Manning's  n and  rock  size  k ; 

A /A  > 0.6 
r 


concrete  surface  yielding  n = 0.01  may  be  regarded  as  a surface 
covered  with  asperities  of  grain  size  equal  to  0.001  ft. 

Method  of  establishing 

the  effective  n of  ejecta 

and  fallback  of  nuclear  craters 

9-  When  the  asperities  of  a channel  are  uniform,  the  Manning's  n 
of  the  channel  is  readily  obtained.  In  the  case  of  ejecta  and  fallback 
of  crater  channels,  the  rock  sizes  are  highly  variable  and  the  deter- 
mination of  n is  not  so  simple.  In  his  experiments  with  pyramidal 

7 

asperities  affixed  to  the  sides  of  a square  closed  conduit,  Fage 
found  that  the  shear  stress  of  the  walls  is  related  to  the  resistive 
force  of  the  pyramids.  He  determined  the  resistive  force  on  the  pyra- 

Q 

mids  by  measuring  the  pressure  on  the  faces  of  the  pyramids.  Sadron 
conceived  that  if  U,  is  the  velocity  at  the  apices  of  the  asperities 

iC 

and  a^  is  the  projected  area  of  the  asperities  in  a plane  normal  to 

the  direction  of  flow,  the  resistive  force  F is 

’ r 


11 


(1.11) 


r 


i 1 


Li 


u,. 


F = Ca  p -f 
r r d 


where  r 

C = the  drag  coefficient 
p = density  of  water 

Next,  Schlichting  found  that  C is  a function  of  the  form  of  the  as- 
perities and  the  concentration  A^/A  . In  principle,  then,  it  would  be 
possible  to  determine  the  shear  stress  of  a surface  covered  with  rocks 
of  varying  sizes  and  concentrations.  As  the  analysis  would  be  very  much 
involved  and  tedious  we  were  not  very  sure  about  the  success  of  this  ap- 
proach. The  method  described  below  has  the  advantage  of  being  simpler 
and  will  be  used  in  the  estimate  of  the  effective  n of  nuclear  crater 
channels. 

Effective  n of  normalized  channels 


10.  Rocks  of  the  surface  of  a crater  channel  are  conceptually 
sorted  according  to  their  sizes  and  are  rearranged  on  the  same  surface 
with  the  provision  that  the  surface  width  everywhere  is  the  same.  This 
is  the  normalized  channel  and  the  arrangement  is  shown  in  Figure  5.  The 
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Figure  5*  Sorted  rocks  of  a crater  channel 


lengths  of  the  rectangular  areas  occupied  by  rocks  of  sizes  k^  k2, 
k^....^  are  1^  12,  l^....,^  , respectively.  Let  the  total  length 


12 


be  L . Manning's  n corresponding  to  the  several  rocks  in  the  separate 
areas  is  n^,  n2,  n^,....,nn  , respectively;  these  values  are  determined 
by  noting  the  rock  sizes.  In  effect  we  now  have  a channel  of  width  B 
and  of  length  L . Imagine  that  a discharge  Q is  related  to  the  depth 
of  waters  H^,  H^,  H^,....,  Hn  . Hq  and  Hn  are  values  at  the 

ends,  and  , say,  is  the  depth  at  the  separation  between  the  rocks 
Assume  that  water  is  large  and  hydraulic  radius  is 
nearly  constant  throughout.  Applying  Manning's  formula 

HQ  - H1  = l^N  N = U2/2.l8RU/3 


k^  and  k2 


H1  - H2  = l2n2N 


H , - H = 1 n N 
n-1  n n n 


H_  - H = ( T l.n. J N 
0 n \i=l  1 V 

Introducing  the  effective  ng  of  the  normalized  channel 


adding 


This  constitutes  the  basis  of  estimating  the  effective  Manning's  n of 
the  crater  channels.  It  is  assumed  that  by  rearranging  the  surface 
rocks  in  the  above  manner,  the  resistance  of  the  channel  thus  formed 
would  be  practically  the  same  as  in  the  original  or  the  initial  channel 
and  that  as  the  rocks  get  smaller  in  size,  their  number  increases. 


Effective  n of  basalt  craters 
and  of  craters  in  glassy  rhyolite 

11.  The  effective  Manning's  n of  craters  in  basalt  and  in 


glassy  rhyolite  was  determined  on  the  basis  of  data  of  size  distribution 
furnished  by  the  engineers  of  U.  S.  Army  Engineer  Nuclear  Cratering 
Group.  This  information  is  contained  in  Figures  6 and  7 and  also  in 
Table  1.  The  procedure  followed  will  be  illustrated  for  the  case  of  the 
basalt  crater.  It  is  imagined  that  the  partitioning  of  the  rocks  is  in 
10  percent  intervals  as  shown  in  Table  2. 

Table  1 

Grain-Size  Distribution  of  Nuclear  Channel  Materials 


Percent  Finer 
by  Weight 


Rhyolite 


Rhyolite 


Weight 

Basalt 

d,  ft 

Adjuste 

100 

6.0 

5.0 

5.0 

90 

3.7 

1.5 

2.0 

90 

3.0 

0.55 

1.0 

70 

2.3 

0.27 

0.50 

6o 

1.8 

0.13 

0.30 

50 

1.1+ 

0.070 

0.19 

4o 

0.90 

0.038 

0.12 

30 

0.58 

0.022 

0.080 

20 

0.30 

— 

0.050 

10 

0.10 

— 

0.036 

0 

” ’ ” 

Parts 

(1) 


(7,  8,  9) 


Steps 

d,  ft 

d3 

N. 

1 

N.d2 

1 

VL 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

100  to  90 

1+.20 

74.0 

1 

17.6 

0.055 

90  to  80 

3.20 

32.7 

2.26 

23.2 

0.073 

80  to  70 

2.60 

17.1 

4.3 

29.O 

0.091 

70  to  60 

2.00 

8.0 

9.2 

37.0 

0.117 

60  to  50 

1.60 

1+.09 

18.1 

46.3 

0.181 

50  to  1+0 

1.10 

1.37 

54.0 

65.2 

0.205 

1+0  to  30 

0.71+ 

0.1+0 

185 

100.0 

0.315 

30  to  20 

0.1+3 

0.08 

920 

20  to  10 

0.20 

0.008 

9200 

318.0 

1.000 

WEIGHT 


The  median  diameters  of  the  steps  are  shown  in  column  3.  Assuming  that 

the  rocks  are  of  like  density,  the  number  of  rocks  of  a given  size  is 

3 

inversely  proportional  to  d . Column  5 denotes  the  relative  number 

of  rocks  of  different  sizes  as  N.  , which  is  inversely  proportional  to 

d . Column  6 gives  the  relative  values  of  the  areas  covered  by  the 

rock  when  placed  close  to  each  other.  In  an  actual  situation,  the  rocks 

of  the  smaller  sizes  would  form  the  base  layer.  It  will  be  assumed  that 

2 

these  rocks  fall  below  30  percent  by  weight.  The  addition  of  ITd 
values  from  parts  1 to  7 inclusive  gives  the  sum  of  318.  Thus,  the 
ratio  N^d  /3l8  gives  the  quantity  1^/L  and  is  shown  in  the  last 
column.  Having  thus  arrived  at  the  discrete  values  of  1^/L  to  be  as- 
sociated with  a given  rock  size  k , we  may  next  proceed  to  obtain  the 
effective  Manning's  n of  the  normalized  crater  using  Equation  1.12. 
This  is  completed  in  Table  3. 


Table  3 


Parts 

k,  ft 

VL 

n 

l.n2/L 

1 

(1) 

(2) 

(3) 

(4) 

(5) 

1 

4.20 

0.055 

0.0453 

1.101  X 10-1* 

2 

3.20 

0.073 

0.0434 

1.482 

3 

2.60 

0. 04l 

0.0420 

1.605 

4 

2.00 

0.117 

o.o4oo 

1.869 

5 

1.60 

0.l8l 

0.0385 

2.380 

6 

1.10 

0.205 

0.0360 

2.660 

7 

0.74 

0.315 

0.0340 

3.660 

14.757  x 10'k 

Application  of  Equation  1.12  yields: 

2 4 

n = 14.76  x io 
e 


n = 0.0384 
e 

the  effective  n of  the  basalt  craters  according  to  this  approach.  It 
appears  that  with  glassy  rhyolite  the  granulometric  distribution  values 
need  to  be  adjusted  prior  to  the  computation.  It  is  inferred  from 
Table  1 that  the  particles  of  the  range  30  percent  finer  by  weight  are 
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too  small.  Discarding  these,  the  new  adjusted  distribution  values  would 
be  as  shown  in  the  last  column  of  Table  1.  Repeating  the  calculations  as 
for  basalt  and  using  the  adjusted  distribution,  the  effective  resistance 
coefficient  for  glassy  rhyolite  is: 


n = 0.0275 
e 


12.  The  effective  Manning  coefficients  for  the  surface  of  nuclear - 

cratered  channels  derived  by  the  method  described  above  are  the  minimum 

ones.  In  the  actual  channels,  the  larger  rocks  are  not  closely  packed 

but  are  randomly  spaced,  and  the  effective  kg  would  be  greater  than 

the  k of  the  closely  packed  rocks.  Again,  whether  the  rocks  are 
s 

closely  packed  or  are  randomly  distributed,  the  rocks  of  smaller  sizes 
would  be  sheltered  by  the  larger  ones.  The  consequence  of  the  latter 
would  be  a reduction  of  the  areas  of  smaller  rocks  exposed  to  the  action 
of  the  flowing  water.  As  a consequence  of  these  two  conditions,  the  ef- 
fective Manning  coefficient  of  an  actual  cratered  channel  surface  would 
be  a factor  larger  than  the  minimum  value  computed  by  the  above  method. 

A modified  form  of  Equation  1.12  is: 


2 

n 

e 


2 

m 


2 

n. 


1. 

i 

L 


(1.13) 


where  m is  a constant  close  to  unity. 

13.  Since  the  theoretical  evaluation  of  m is  not  practical  it 

is  necessary  to  resort  to  experimental  evidence.  In  recent  experiments 

o 

on  the  flow  of  water  over  graded  beds,  three  sizes  of  gravel  were  em- 
ployed with  gradations  as  shown  in  Table  4.  In  forming  the  bed,  first 
the  larger  rocks  were  randomly  placed  and  successively  the  small  sizes 
were  added.  At  the  ends  the  beds  were  of  two  layers.  We  have  theoret- 
ically determined  the  effective  n^  for  these  gravel  surfaces  by  the 
above  methods  and  also  from  the  flow  measurements.  The  steps  in  the 
determination  of  ng  for  the  gravels  of  gradation  1 are  shown  in 
Tables  5 and  6.  Here,  also,  the  gravels  30  percent  finer  by  weight  are 
assumed  to  be  in  the  underlayer.  The  theoretically  computed  coefficient 


Table  4 

Gradation  of  Gravels  in  the  WES  Experiments' 


Percent 

Finer  by 

Gradation,  d , ft 

Weight 

1 

2 

3 

100 

0.124 

0.0624 

0.0312 

90 

0.101 

0.0525 

0.0231 

8o 

o.oo4 

0.044 

0.0175 

70 

0.071 

0.037 

0. 01 4l 

6o 

0.06l 

0.031 

0.0122 

50 

0.048 

0.026 

0.0103 

40 

0. 040 

0.022 

0.0082 

30 

0.032 

0.019 

0.0075 

20 

0.027 

0.016 

0.0055 

10 

0.022 

0.011 

0.0043 

0 

0.019 

0.008 

0.0039 

Table  5 


Percent 

Finer 

by 

Weight 

d,  ft 

d3 

N. 

1 

d2 

N.d2 

1 

h/L 

100-90 

11.25  x 10-2 

1420  x 10-6 

1.00 

126  X 10“U 

126  x 10-U 

0.079 

90-80 

9-25 

791 

1.79 

85.5 

153 

0.092 

80-70 

7-75 

465 

3.05 

60.0 

183 

0.109 

70-60 

6.20 

238 

5.96 

38.4 

228 

0.136 

60-50 

5.45 

162 

8.76 

29.7 

260 

0.155 

50-40 

4.40 

85 

16.70 

19.4 

325 

0.194 

40-30 

3.60 

46.6 

30.5 

13.0 

396 

0.237 

30-20 

2.95 

25.7 

55.2 

8.7 

462 

20-10 

2.45 

14.7 

96.8 

6.0 

580 

10-0 

2.10 

9.26 

153 

48 

730 

Table  6 


Parts 

k,  ft 

VL 

n 

2 

n 

n2li/L 

1 

0.1125 

0.077 

0.0250 

6.26  X 10~k 

0.482  x io-4 

2 

0.0925 

0.092 

0.0240 

5.78 

0.532 

3 

0.0775 

0.109 

0. 0233 

5.44 

0.594 

4 

0.0620 

0.136 

0.0225 

5.05 

0.688 

5 

0.0545 

0.155 

0.0220 

4.85 

0.752 

6 

0.0440 

0.194 

0.0215 

4.62 

O.890 

7 

0.0360 

0.237 

0.0205 

4.20 

1.00C 

n2  = 4.946 

-4 

x 10  ; 

n = 0.0222 

: _4 

4.,.u  x 10 

--  . «. 


ne  equals  0.0222  for  the  gravel  surface  of  gradation  1 and  0.0199  and 

0. 0171  for  gradations  2 and  3,  respectively.  In  Reference  9 the  values 
of  the  ratio  U/U#  and  the  depths  for  each  test  are  recorded.  Using 
these  values,  we  have  determined  hy  Equation  1.7  the  n for  every  test 
run.  There  is  considerable  disparity  in  the  values  of  n from  the  tests 
with  a given  distribution  of  rocks,  but  by  averaging  the  observed  values, 
these  are  ng  = 0.0280,  0.02^3,  and  0.0213  for  surfaces  of  gradations 

1,  2,  and  3,  respectively.  Comparing  these  with  the  theoretical  values 
shown  and  using  m for  the  ratio  of  the  observed  to  the  computed,  we 
have  m = 1.25,  1.22,  and  1.2b  for  the  gravels  of  gradations  1,  2,  and 
3,  respectively.  Hence,  Equation  1.13  can  be  modified  as: 


Applying  this,  a better  value  of  resistance  coefficient  for  the  basalt 
craters  is 

n = 1.2U(0.038U)  = O.OU76 
e 

and  for  glassy  rhyolite: 


n = 1.2^(0.0275)  = 0.03^2 
e 


lU.  An  alternate  procedure  in  the  last  problem  would  be  to  con- 
sider the  effective  k of  a given  gradation  of  rc  ks.  In  the  WES  test, 
9 8 

Brown  has  determined  the  kg  for  each  run.  There  are  large  disparities 
between  the  individual  values  from  the  tests  and  a bed  of  a given  gra- 
dation of  rocks.  If  averages  are  taken,  k is  equal  to  0.233,  0.0968, 

s 

and  0.0613  ft.  One  should  inquire  as  to  the  ratio  of  k in  the  aver- 

s 

age  value  of  rock  sizes  for  a given  gradation.  The  average  value  of 
the  rock  sizes  d may  be  determined  as 


or 


_ dl. 

d « Z (1.15) 


as  it  is  best  to  do  the  averaging  on  the  basis  of  areas  covered  by  the 


rocks.  One  finds  the  mean  values  equal  to  0.0598,  0.0310,  and  0.0125  ft 
for  rocks  of  the  gradations  1,  2,  and  3.  Thus 


As  an  average,  then. 


— = 3.89  for  gradation  1 
d 

= 3.13  for  gradation  2 
= li.90  for  gradation  3 


kg  = 3. 97  d 


(1.16) 


Resorting  to  the  values  of  d and  1^/L  shown  in  Table  2,  d = 1.687  ft 
for  the  basalt  rocks.  Hence,  using  Equation  l.l6,  kg  = 6.69  ft 
and  turning  to  Equation  1.8,  the  channel  effective  n is 


ne  = 0.0M5 


which  agrees  fairly  well  with  the  calculation  by  Equation  l.lU.  In  the 
same  manner,  one  finds  for  the  rhyolite  channel 


n = 0.032J+ 
e 


Loss  of  Energy  in  an  Isostatic  Area 


Losses  in  isostatic  flow  areas 

15.  The  term  isostatic  refers  to  areas  of  the  flow  where  a mean 
movement  out  of  the  area  is  absent,  although  there  are  local  movements 
from  point  to  point.  In  these  areas,  pressures  are  uniform  and  the 
movements  are  essentially  cyclic  induced  by  the  tangential  stresses  of 
the  main  flow  turbulence.  This  situation  arises  when  a current  is 
deflected  or  a current  moving  in  a constant  direction  experiences  a 
sudden  expansion.  In  such  areas  there  is  loss  of  energy  in  the  main 
current,  due  partly  to  the  maintenance  of  the  eddying  motion  and  to 
the  turbulent  dissipation  in  the  waters  adjacent  to  the  isostatic  area. 


V *ri'-  f' 


Energy  loss  in  Tollmien  free  jet 

16.  When  a current  meets  and  runs  over  a body  of  still  water,  a 
zone  of  turbulent  waters  is  created  as  shown  in  Figure  8.  The  zone  is 
limited  between  two  straight  lines  oy1  and  oy2  showing  that  the 
width  of  the  turbulent  zone  increases  linearly  with  distance  measured 
from  the  point  where  the  current  meets  the  still  water.  At  any  point  on 


Figure  8.  Free  turbulence  between  current  and 
still  water  (semijet) 

the  line  oy^  , the  longitudinal  velocity  equals  the  initial  velocity 

u^  , of  the  current.  At  any  point  on  the  line  oy2  , the  longitudinal 

velocity  vanishes  but  there  is  a lateral  movement.  The  line  oy^  is 

the  zero  streamline  with  an  angle  of  inclination  a with  the  horizontal. 

The  water  below  oy  is  in  the  isostatic  area.  The  theory  of  this  type 
0 10 

of  mixing  has  been  developed  by  Tollmien  and  it  can  be  shown  that 
tan  a = 0.19  . A single  streamline  in  this  area  closes  at  infinity  and 
the  flow  below  oy^  is  virtually  cyclic.  Energy  lost  from  the  main 
current  is  expended  in  maintaining  the  cyclic  motion  below  oy^  and  in 
the  turbulence  of  the  zone  above  oy 2 . Analysis  shows  that  the  jet 
dissipation  (loss  of  energy  per  unit  time)  amounts  to 
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; 

i 

i 

i 

i 
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\ 

c 

i 

I 


El  = 0.0468  p 2^Ax  (1.17) 

where  is  the  area  of  contact  between  the  isostatic  zone  of  cyclic 

motion  and  the  main  current.  The  losses  encountered  in  miter  bends  and 
in  sudden  expansions  are  nearly  in  accordance  with  Equation  1.17  as 
shown  by  Sager  and  others.1'1'  Thus  this  should  prove  to  be  of  value 
when  considering  the  additional  losses  in  cratered  channels  owing  to 
the  sudden  changes  of  cross  sections. 

Formulation  of  Losses  of  Nuclear-Cratered  Channels 

17.  A channel  created  by  nuclear  devices,  placed  equidistantly 
along  a practically  straight  line,  will  consist  of  a succession  of  chan- 
nel segments,  an  individual  segment  having  the  planar  form  as  shown  in 
Figure  9*  The  segment  in  one  part  has  divergent  boundaries  and  in  the 
remaining  part,  convergent  boundaries.  A segment  of  this  configuration 


R = AVERAGE  HYDRAULIC  RADIUS 
X = AVERAGE  WETTED  PERIMETER 
U = AVERAGE  CURRENT  VELOCITY 


Figure  9-  Notation  diagram  of  cratered  segments 

will  be  referred  to  as  a cratered  segment.  Despite  the  fact  that  the 
devices  creating  the  channel  are  arranged  equidistantly,  the  resulting 
crater  segments  may  not  always  have  equal  divergent  and  the  convergent 
parts.  Let  L be  the  total  length  of  the  cratered  segment  and  be 

the  length  of  the  divergent  part.  If  the  flow  were  reversed,  the  length 
of  divergent  part  would  be  L - . Taking  the  point  of  view  adapted  in 

22 
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the  previous  sections,  it  is  assumed  that  a flow  separation  occurs  in  the 

l 

divergent  part  of  the  cratered  segment.  This  is  part  of  an  isostatic 
area,  and  in  this  area,  the  loss  of  energy  per  unit  time  is: 


\ = C | U3A1 


(1.18) 


where 

= jet  dissipation  due  to  free  turbulence  in  a nuclear  segment 

C = an  experimentally  determined  constant  depending  on  the  form 
of  the  segment 

U = the  mean  velocity  averaged  along  the  segment  length  L 
= the  total  area  of  the  separation  surface 

In  additicn,  there  is  the  rate  of  loss  in  the  segment  induced  by  the 
solid  boundaries 


1 ■? 

e2  = — PUJ(A  - A1) 


(1.19) 


where 


- the  coefficient  of  friction  due  to  the  boundary  asperities 
A = the  wetted  wall  surface 

If  n1  is  Manning's  n for  the  surface,  from  Equation  1.7 


, 29.16  2 
*i " 7173  “1 


(1.20) 


where  R is  the  mean  hydraulic  radius  of  the  cratered  segment. 

l8.  Let  Q be  the  discharge  and  AH  the  fall  of  the  energy  line 
along  the  segment.  Hence, 


pgAHQ  = E1  + E2 


gAHpQ  = X1  | U3(A  - A±)  + C | U3A1 


(1.21) 


Introducing  the  average  wetted  perimeter  x 


Q = xRU,  A = Lx,  A1  = ^x 


Substituting  these  in  the  above  and  dividing  by  pQL  one  has 


(1.22) 


This  may  be  written,  introducing  X 2 , the  effective  resistance  coeffi- 
cient of  the  cratered  segment,  as 


g 


AH 

L 


X ^ 
2 2R 


where 


One  may  also  write  this  as 


where 


(1.23) 


(1.21) 


(1.25) 


(1.26) 


19.  Equation  1.25  above  refers  to  a single  cratered  segment.  In 
the  derivation,  use  was  made  of  the  average  value  of  the  mean  currents 
in  various  channel  cross  sections.  This  is  permissible  if  it  is  assumed 
that  the  variation  of  the  cross  sections  and  of  the  hydraulic  radii  be- 
tween extremes  is  linear  with  distance,  flow,  if  it  is  further  assumed 
that  the  quantities  Jl^/L  and  C are  practically  constant  over  the 
channel  length.  Equation  1.21  may  be  assumed  to  apply  to  a long  cratered 
channel  consisting  of  many  cratered  segments,  and  AH/L  would  denote  the 
gradient  of  energy  line  of  the  long  channel.  The  interpretation  of 


Equation  1.23  is  that  once  the  average  quantities  l^/L  and  N are 
known,  the  effective  X 2 or  the  Manning's  n 2 of  the  channel  can  be 
determined.  For  determination  of  these  quantities  in  special  cases, 
reference  will  be  made  to  the  WES  channel  experiments."^ 

Results  from  WES  Channel  Experiments 

20.  The  nature  of  the  increased  losses  due  to  the  divergence  and 
convergence  in  a cratered  channel  was  examined  in  three  channels,  all 
having  the  same  length,  the  same  average  cross  section  of  1.206  sq  ft, 
and  the  same  average  hydraulic  radius  of  0.0378  ft.  The  surface  asperi- 
ties of  all  the  channels  were  the  same.  One  of  the  channels  was  of  uni- 
form cross  section  throughout.  The  other  two  channels  consisted  of 
cratered  segments  with  form  specifications  indicated  by  the  coordinator 
of  the  WES  channel  experiments.  Major  Robert  G.  Bening,  U.  S.  Army 
Engineer  Nuclear  Cratering  Group,  Livermore,  California.  The  experi- 
ments showed  that  the  losses  of  the  two  cratered  channels  were  practi- 
cally of  like  value. 


—■  = 4.65  x 10  ^Q2 

and  for  the  channel  of  uniform  cross  section 

^ = 1.36  x 10~^Q2 

Using  Manning's  formula  for  both  channels 

2 2 
AH  _ n U 

L 2.22  TT/3 


for  the  cratered  channel  (1.23  bis) 


and 


AH 


'2  U 


g L 2 R 


and 
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(1.20  bis) 


l/o  p 

X^R  ' = 29.16  for  the  smooth  channel 


one  finds  that  for  the  uniform  channel 


= 48.1  x 10 


-4 


and 


n±  = 0.0108 


For  the  cratered  channels 


X2  = 164.5  x 10 


-4 


and 


n2  = 0.0201 

Accordingly,  the  expansions  and  the  contractions  do  increase  the 
Manning's  n to  twice  the  value  of  the  uniform  channel. 

21.  The  area  of  the  separation  of  the  flow  in  the  cratered  chan- 
nel was  examined  visually  through  the  introduction  of  dye  in  one  of  the 
cratered  channels  for  the  two  directions  of  flow,  regular  and  reversed, 
and  for  two  discharges.  Only  the  surface  dimensions  of  the  separations 
could  be  noted  with  certainty.  Let  g-| m denote  the  length  of  the 
separated  area  in  the  direction  of  flow,  A denote  the  maximum  width 
of  this  same  area,  and  denote  the  maximum  width  of  the  segment  (see 

Figure  10).  The  observed  values  of  /L  show  some  variations:  for 
one  discharge  the  average  values  were  O.698  and  O.678  with  reversed  flow; 
for  another  discharge  the  corresponding  values  were  0.638  and  0.687. 

Also,  there  were  small  variations  in  the  values  of  2A/B  from  segment 

m 

to  segment:  for  one  discharge  the  average  values  were  0.263  and  0.288; 
for  the  other  discharge  the  corresponding  values  were  0.285  and  0.354. 
These  show  the  extent  of  the  separations.  The  same  dye  technique  could 
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Figure  10.  Notation  diagram  of  surface  separations 
in  cratered  segments 

not  be  extended  to  points  below  the  surface.  The  impressions  gathered 
indicated  that  the  separation  became  smaller  with  increased  depth.  As 
the  extent  of  the  separations  should  depend  on  the  inclinations  of  the 
boundaries,  the  depth  dependence  should  be  noted  for  the  channels  of  the 
investigation.  The  inclination  6 varies  with  depth.  Values  for  the 
different  crater  segments  do  vary  somewhat.  Without  making  a distinc- 
tion between  the  inclinations  in  the  divergent  and  convergent  parts, 

the  mean  of  0 at  the  surface  is  9*95  degrees  and  0^  at  the  bottom 
s b 

6.62  degrees  for  the  two  channels. 

Establishment  of  Cratered  Channel  Formula 

22.  Determinations  of  the  quantities  £^/L  and  N of  Equa- 
tion 1.25  will  be  made  on  the  basis  of  the  results  obtained  in  the  WES 
cratered  channel  experiments.  The  quantity  ^ can  be  related  to  the 
length  & of  the  separation  area  observed  at  the  surface  in  the  fol- 
lowing manner.  One  may  consider  the  wetted  perimeter  of  a cross  section 
in  the  general  area  of  the  separation  (Figure  10).  Let  s be  the  arc 
length  measured  from  the  lowest  point  in  the  section  along  the  rigid 
boundary.  Place  2s^  = x • Let  £ be  the  length  of  the  separated  area 
at  the  level  of  s . Thus,  , the  total  area  of  the  separating 

surface,  is 
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As  mentioned  previously,  the  length  of  separation  decreases  for  the 
deeper  points.  Nov,  assume  that 


Accordingly,  since  2s^  = x 


Al  = 


1 + 


m ^lm 


Recalling  the  definition. 


, it  is  seen  that 


1 + m 


and 


L 1 + m L 


23.  The  experimental  results  (shown  in  Tables  4 and  5)  show  that 
the  average  value  of  Jl^/L  from  four  sets  of  separate  observations  is 
0.668.  If  one  places  coefficient  m = 1.5  , then 


^1  _ 0.668 

L “ 2.5 


0.268 


Equation  1.25  then  becomes 


\2  = 0. 732X1  + N 

In  the  WES  channel  tests  with  the  channel  of  uniform  cross  section 

X1  = 48.1  x io~h 
and  in  the  cratered  channels 
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X2  = 16U.5  X 10  11 

With  these  values  it  is  found  that 

N = 129.3  x 10-1* 


with  m = 3.0 


1 

= 0.223 

\2  = 0.77TA1  + N 

and 

N = 127.3  x 10"1* 

Again,  if  m = 1.0  , 


E-  0-33l. 

Ag  = 0.666A1  + N 


and 


N = 132.3  x 10' 


-It 


(1.27) 


(1.28) 


(1.29) 


2l|.  In  the  above  there  are  three  expressions  to  evaluate  the 
effective  resistance  coefficients  of  a cratered  channel.  It  is  diffi- 
cult to  make  a selection.  Before  a selection  is  made,  however,  it  is 
necessary  to  ascertain  if  these  expressions  could  be  valid,  at  least 
approximately.  The  criterion  is  that  the  losses  thereby  indicated  are 
of  the  order  met  within  the  turbulent  isostatic  zone  of  a semijet.  It 
will  be  recalled  that  the  latter  loss  is  given  by 


C = 0.01*68 


(1.17  bis) 


29 


With  the  case  of  the  cratered  segments,  correspondingly,  the  expression 
is 

E,  = c | u\  (1.18  bis) 

and 


c = N 


(1.26  bis) 


Introducing  values  of  N and  L/i^  from  Equations  1.27,  1.28,  and 
1.29  gives 

m = 3,  C = 0.0572 


m = 3/2,  C = 0.01*83 


m = 1,  C = 0.0393 

which  show,  indeed,  that  the  losses  thus  calculated  for  a cratered 
segment  are  of  the  same  order  of  magnitude  as  for  a semijet.  We  shall 
adopt  the  expression  of  loss  corresponding  to  m = 3/2  . Accordingly, 
the  expression 

A2  = 0.732A1  + 129.3  x io_U  (1.30) 

is  adopted  for  the  effective  coefficient  of  friction  of  a cratered 
channel  segment  of  the  general  form  tested  in  the  model  channels  at  WES. 


Cratered  Channel  Loss  Formula  in  Terms  of  Manning's  n 

25.  The  cratered  channel  loss  expression 

X2  = 0.732A1  + 0.0129  (1.31) 

involving  the  roughness  coefficients  may  be  changed  into  one  involving 
the  Manning's  coefficients  n2  and  n^  . Using  the  form  in  Equation 
1.20,  measuring  R in  feet. 


(1.32) 


P ll  P ]l  1 / Q 

x 10  = 0. 732n^  x 10  + 1*.1*3R'  ' 


where 

n2  = the  effective  Manning's  n of  the  cratered  channel 

n^  = the  Manning's  n if  the  channel  did  not  have  expansions  and 
contractions 

R = the  average  value  of  the  hydraulic  radii 


The  above  expression  indicates  that  the  increase  of  the  effective  Man- 
ning's n of  a cratered  channel  over  that  of  a uniform  channel  would  be 
more  pronounced  with  smaller  surface  asperities  and  with  greater  hydrau- 
lic radii. 

26.  It  is  envisaged  that  the  formula  given  could  not  be  applied 
to  all  cases  of  cratered  channels  but  only  to  those  having  form  char- 
acteristics similar  to  the  channels  studied  at  WES.  As  a criterion  of 
form,  one  may  take  the  angle  of  divergence  0 of  the  channel  lines  at 

the  water  surface  and  those  at  the  bed  bottom.  Let  0 refer  to  the 

s 

surface  expansions  or  contractions  and  0^  to  the  bottom  expansions  or 
contractions.  Taking  average  values  of  these  quantities  expressed  as 
radians  from  all  the  cratered  segments  in  channel  2 and  in  channel  3: 

0g  = 9.95  degrees  and  0^  = 6.62  degrees  . Note  than  0^  is  smaller 
than  0g  . This  may  be  the  reason  that  the  separation  areas  in  the  ex- 
panding sections  are  greater  and  more  significant  at  the  water  surface. 

27.  An  alternate  approach  to  the  problem  is  one  proposed  by 
Mr.  J.  H.  Douma  of  the  Office,  Chief  of  Engineers,  U.  S.  Army.  The 
losses  are  visualized  as  arising  from  expansion.  The  idea  was  followed 
previously  by  Sager  and  others’'-'*'  with  the  result  that 


n2  x loh  = x id*  + 3-99R1/3  (1.33) 

One  may  call  this  the  Douma  formula.  The  difference  between  the  main 
formula  (Equation  1.32)  and  the  Douma  formula  (Equation  1.33)  is  small. 
Indeed,  theoretically.  Equation  1.33  represents  the  limiting  form  of 
Equation  1.32.  The  Douma  formula  could  have  been  obtained  directly  from 
Equation  1.2h  if  one  had  supposed  that  in  this  equation  H^/L  vanishes 
and  C becomes  infinite.  Interpreting  physically,  one  would  have 
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assumed  that  the  separation  areas  responsible  for  added  losses  are  in- 
finitely small  and  the  dissipation  rate  infinite.  This  is  a difficulty 
avoided  in  the  derivation  of  Equation  1.32,  which  is  suitable  for  the 
estimate  of  the  effective  resistance  of  cratered  channels. 


Effective  Resistance  of  Nuclear-Cratered  Section 


28.  In  August  1968  the  Nuclear  Cratering  Group,  Corps  of  Engineers, 
furnished  data  on  dimensions  of  a long  section  of  a hypothetical  nuclear- 
cratered  channel.  These  data  will  be  utilized  to  illustrate  the  manner 
of  determining  the  effective  resistance  using  Equation  1.32.  The  size 
distribution  is  given  in  Table  7,  the  material  of  the  section  being 
basalt,  rhyolite,  or  clay  shale. 


Table  7 

Nuclear  Excavation  Data 


Explosive 

Location 

Explosive 

Spacing 

ft 

Channel 

Depth 

ft 

Channel 

Width 

ft 

Channel  Area 
in  Cross  Section 
sq  ft 

Hydraulic 

Radius 

ft 

Channel 

Length 

ft 

1 

565 

336 

1193 

227,928 

165 

0 

2 

652 

1*18 

1436 

336,969 

201 

565 

3 

739 

312 

1188 

216,573 

160 

1,217 

1* 

652 

386 

1407 

312,625 

193 

1,956 

5 

652 

374 

1307 

222,002 

182 

2,608 

6 

823 

L6l 

1627 

426,417 

226 

3,260 

7 

907 

330 

1300 

253,142 

172 

4,083 

8 

823 

361* 

1402 

299,075 

187 

4,990 

9 

739 

420 

1507 

S62.169 

208 

5,813 

10 

739 

3l*0 

1271 

251,012 

172 

6,562 

11 

823 

41*0 

1566 

' 392,895 

217 

7,291 

12 

823 

293 

1188 

207,155 

156 

8,114 

13 

652 

4 61 

1627 

426,417 

226 

8,937 

14 

565 

374 

1307 

276,664 

182 

9,589 

1 

335 

1193 

227,928 

165 

10,154 

29.  Let  L be  the  length  of  a segment,  R and  A the  average 
hydraulic  radius  and  the  average  cross  section,  respectively,  of  the 
segment,  Q the  discharge,  and  AH  the  fall  of  the  water  surface  be- 
tween the  ends  of  the  segment.  Now 


1 


AH  = n NL 
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N = 0.1*58  — & 


2 

2^73 


where  n is  the  Manning's  coefficient  of  the  segment.  Dropping  the 
subscript  2,  Equation  1.32  becomes 


n2  x IQ1*  = 0.732n2  * IQ**  + 1*.1*3R1/3 


The  fall  of  the  water  surface  in  the  individual  segments,  in  the  order 
of  the  segments,  is 


H 

H 


1 

2 


H - 

n 


n+1 


Lin?»i 


^J2n2N2 


L n2N 
n n n 


Similarly,  introducing  the  effective  Manning's  n for  the  entire  section 
of  length  Lt  , 

Hn+1  - “l  ’ " 


where  N is  computed  as 


N - 0.1*58  =5^73 


where  A and  R are  the  average  cross  section  and  the  average  hydraulic 
radius,  respectively,  of  the  entire  channel.  Since  the  energy  loss  in 
the  entire  channel  equals  the  sum  of  the  losses  in  the  segments: 


and  hence 
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2 

n 

e 


(1.34) 


If  and  are  the  hydraulic  radii  at  the  ends  of  a segment  and 

R2  the  hydraulic  radius  at  the  cross  section  of  the  largest  area,  and 
and  L2  are  the  lengths  of  converging  and  diverging  parts,  then  the 
average  hydraulic  radius  of  a segment  i is 


f r2  + R3' 


Z 

L 


Similarly  the  average  cross  section  of  a segment  i is 


The  average  hydraulic  radius  R and  the  average  cross  section  A of 
the  entire  section  are  evaluated  as: 

n L. 

R = V R. 

A 1 lt 


and 


n 

A=  I 


A. 

l 


L. 

l 


i=l  * T 

Equation  1.3*+  requires  that  for  the  evaluation  of  the  effective  n 

e 

the  entire  channel  the  quantity 


of 


relating  to  a segment  of  position  i be  determined.  For  the  hypotheti- 
cal channel  of  the  dimensions  shown  in  Table  7,  the  steps  of  the  compu- 
tations are  shown  in  Table  8. 

30.  To  complete  the  computation  of  ng  through  Equation  1.34, 
the  n2  for  each  segment  must  also  be  known.  As  will  be  recalled, 
this  is  done  using  Equation  1.32.  For  a basalt  hypothetical  channel, 
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Table  8 


Steps  in  the  Evaluation  of  n 


Segment 


Segment 


Average 


Segment 


Segment 


n_L  = 0.0476  . This  would  be  the  appropriate  Manning's  n had  the 

channel  been  of  uniform  cross  section  throughout.  Then  the  appropr: 

values  n„  on  the  basis  of  Equation  1.32  are  given  in  the  fourth  c< 
2 _ 2 
of  Table  9.  The  sum  of  the  entries  in  the  last  column  give  ng  x 1 

thus  the  effective  Manning's  n of  the  hypothetical  channel  of  the 


I 


Table  9 

Coefficient  of  Basalt  Section 


R1/3 

1 

4.43rV3 

2 1* 

n*  x 10 

u . 

1 

lt 

(a)  • (b) 

5.63 

2k.  9k 

1*7.69 

0.1523 

7.263 

5-67 

25-12 

1+7.87 

0.1874 

8.970 

5.86 

25.96 

1*8.21 

0.1109 

5.335 

5.75 

25-47 

48.22 

0.2346 

11.313 

5.75 

25.1*7 

48.22 

0.1463 

7.044 

5.75 

25.1*7 

48.22 

0.1892 

9.123 

Sum 

58.02 

= 58.02 

n = 0.0762 
e 

2 4 

n x 10 
e 


dimensions  shown  in  Table  7 is  ne  = 0.0762  . Making  similar  computa- 
tions for  channels  of  rhyolite  or  clay  shale,  the  results  are  shown 
below. 

n,  n 


Material 

Basalt 
Rhyolite 
Clay  shale 


0.0476 

0.0332 

0.010 


0.0762 
0. 0609 
0.  0514 


An  obvious  deduction  from  this  summary  is  that  divergence  and  convergence 
in  segments  considerably  increase  the  resistance  factor  of  the  channnels. 

Estimate  of  Roughness  of  Channels  with  a Muck  Bottom 


31.  According  to  Mr.  Eden,  formerly  of  SAJ,  the  distribution  of 
grain  size  of  the  Atlantic  muck  is  100  percent  less  than  0.1  mm,  65  per- 
cent less  than  0.01  mm,  and  35  percent  less  than  0.002  mm  (see  Figure  24). 
The  estimate  of  roughness  having  a material  this  fine  raises  some  diffi- 
cult questions.  If  it  is  assumed  that  no  ripples  or  dunes  are  present, 
the  boundary  surface  could  be  hydrodynamically  smooth.  In  order  that 
this  be  the  case,  the  criterion 


(1.35") 


need  be  satisfied.  Here  U-  is  the  shear  velocity  Ji  /p  , t the 

o o 

wall  stress,  p the  density  of  water,  v the  kinematic  viscosity,  and 
d the  grain  diameter.  With  this  condition  satisfied,  the  mean  current 
U prevailing  in  a trapezoidal  channel  would  be 


jf-  = 3.25  + 5.75  log 

U# 


(1.36) 


One  may  approximate  this  by  a power  expression  of  the  form 


(1.37) 


The  coefficient  K , a dimensionless  number,  may  be  determined  by  equat- 
ing the  right-hand  sides  of  Equations  1.36  and  1.37»  for  a selected 
value  of  U#R/v  . The  values  of  K obtained  in  this  manner  are  shown 
in  Table  10.  Comparing  Equation  1.37  with  Manning's  law.  Equation  1.3, 
the  effective  sand  roughness  appropriate  for  hydrodynamical  smooth 
surface  is 


Table  10 


U*R 

UR 

4 

Ki 

K2 

V 

V 

K 

1 x 10!l  ' 

1 x ioJ 

1 x 10 

1 x 10? 

1 x 10^ 

1 x 10 1 

1.1+8  x 
2.05  x 

2.65  x 
3.20  x 
3.75  x 
U. 35  * 

10?  6.85 

10^  6.50 

10;  5.75 

10°  U.70 

loi  3.77 

10  2.97 

0.0382 

0.01+05 

0.01+56 

0.0560 

0.0697 

0.0881+ 

0.0501+ 

0.0532 

0. 0601 
0.0700 
0.081+3 
0.1032 

< 

K 

i. 

(1.38) 

Using  Equation 

1.8 

rrf)"1'6 

(1.39) 

where 
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ilft.lii KVT  »i  » 


0.032U  x 8.12 

K 


(1.39) 


K 


This  establishes  the  Manning's  n for  the  hydrodynamically  smooth 
surfaces.  The  dependence  of  the  coefficient  K on  the  flow  Reynolds 
number  based  on  shear  velocity  is  shown  in  Table  10. 

32.  The  effective  sand  roughness  depends  on  the  channel  roughness 
boundary  stress.  The  determination  of  this  quantity  in  a smooth  channel 
is  somewhat  uncertain  if  the  flow  of  the  channel  is  unsteady.  In  simple 
terms,  the  equation  of  motion  of  this  flow  would  be 

-u2 

3U  3h  U» 

3t  g 3x  R 


Now,  if  the  resistance  term  (the  right-hana  side)  is  large  in  comparison 
with  the  acceleration  term,  the  boundary  stress  is  determined  by  the 
slope  of  the  water  surface.  Thus,  if  L is  the  length  of  the  inter- 
oceanic  channel,  R is  the  hydraulic  radius,  and  H is  the''7 semirange 
of  the  Pacific  tide,  then  the  limiting  value  of  U#  for  steady  flow  is 
« 

U*  = 


This  could  be  true  only  when  the  tidal  period  extends  into  days.  Since 
the  acceleration  term  in  general  is  not  negligible  with  respect  to  the 
resistance  term,  one  may  represent  this  by  writing 


U 


* 


(1.1*0) 


where  m is  a fractional  number.  Taking  L as  50  miles,  R = 50.6  ft 
which  is  the  hydraulic  radius  of  a trapezoidal  channel  with  dimensions 
of  600  ft  surface  width  and  60  ft  depth  with  an  embankment  slope  of  1 on 
1.5,  one  may  evaluate  Manning's  coefficient  n for  a set  of  values  of 
m and  the  energy  loss  H using  Equation  1.39  with  U#  from  Equation 
1.1*0  and  assuming  that  v = 1 x 10  ft  /sec.  The  result  is  the  follow- 
ing set  of  values  for  n . 
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H,  ft 


m 

10 

2.5 

0.62 

i 

0.0129 

0.0133 

0.0136 

1/2 

0.0133 

0.0136 

0.0139 

1/1* 

0.0136 

0.0139 

0.0156 

1/8 

0.0139 

0.0156 

0.0160 

These  results  indicate  that  in  the  part  of  a canal  with  a muck  boundary 
the  effective  Manning  coefficient  varies  slightly  over  a tidal  cycle. 

33.  An  alternate  treatment  of  this  problem  would  be  to  express 
the  effective  n in  terms  of  mean  current  instead  of  mean  shear  ve- 
locity. From  Equation  1.37 

u.  = r6/7  (!f)'1/T  u 

Inserting  this  in  Equation  1.38 

Hr  - (^)6  *6/7  W6n 

and  introducing  in  Equation  1.8 

„ . K2 

where  Kg  = 0.2631K-6/T 
The  dependence  of  K 2 on  the  Reynolds  number  UR/v  is  shown  in 
Table  10.  The  variation  of  with  Reynolds  number  is  appreciable. 

In  a given  channel  the  Manning's  n decreases  when  the  velocity  of 
current  is  increased.-  Considering  again  the  same  channel  as  described 
in  the  preceding  paragraph,  R = 50.6  ft  and  v = 1 x 10  ft  /sec,  the 
following  set  of  n values  is  calculated  for  various  current  velocities. 


ft /sec 

n 

U,  ft/sec 

n 

0.5 

0.0154 

6.0 

0.0136 

1.0 

0.0152 

8.0 

0.0131 

2.0 

0.0146 

10.0 

0.0130 

4.0 

0.0139 

34.  It  is  also  necessary  to  establish  the  limiting  value  of  the 


(1.1*1) 


(1A2) 


(1.1+3) 


39 


current  above  which  the  channel  walls  cease  to  be  hydrodynamically 
smooth.  To  determine  the  limiting  velocity,  the  value  of  U#  from 
Equation  l.Ul  is  introduced  into  Equation  1.35.  Thus,  the  criterion 
for  hydrodynamic  smoothness  of  a channel  is 

7/6 

U.02K  1 

v — 


— < k.OZK^j 


U.kk) 


where  d is  the  grain  size  dg^,  which  for  the  Atlantic  muck  is  0.01  mm 


on  3.28  x io  ft.  Inserting  this  in  the  above,  taking  R = 50.6  ft, 
v = 1 x 10  ft  /sec  and  obtaining  the  value  of  K from  Table  10  as  3, 
the  maximum  value  of  the  current  above  which  the  channel  ceases  to  be 
hydrodynamically  smooth  is  Ul. 5 ft/sec.  Thus,  the  channel  may  be  con- 
sidered to  be  hydrodynamically  smooth  under  all  flow  rates.  On  this 
basis,  the  average  value  of  the  Manning's  n of  a channel  in  Atlantic 
muck  with  no  ripples  or  dunes  during  a tidal  cycle  as  deduced  from  the 
results  of  the  preceding  paragraphs  would  be  about  0.014.  This  also 
applies  to  clay  shale  with  the  grain-size  distribution  shown  in  Figure  11. 


PART  II:  MECHANICS  OF  BED  FEATURES  IN  ALLUVIAL  CHANNELS 


Entrainment  of  Sand 


35.  In  alluvial  channels,  where  the  bed  is  composed  of  loose 
sand,  the  resistance  encountered  for  a given  flow  is  much  higher  than 
what  may  be  inferred  from  the  drag  on  the  individual  grains.  This  in- 
crease has  been  attributed  to  the  form  resistance  of  the  bed  features. 

The  low  Manning  n values  derived  in  previous  paragraphs  for  channels 
in  Atlantic  muck  were  derived  under  the  assumption  that  no  bed  features 
were  present.  In  natural  alluvial  streams  characterized  by  bed  features, 
the  grain  size  is  about  0.1  mm.  It  is  not  clear  if  bed  features  should 
be  expected  with  grain  sizes  as  small  as  in  Atlantic  muck  or  clay  shale 
of  the  same  size;  this  point  should  be  examined.  The  answer  may  pos- 
sibly be  obtained  from  the  complex  hydraulics  of  movable-bed  material. 

The  literature  in  this  field  of  hydraulics  is  extensive.  Nevertheless, 
it  may  suffice  to  focus  attention  only  on  the  form  resistance  of  bed 
features  of  natural  channels,  laying  stress  on  laboratory  studies,  and 
to  develop  proper  formulations  to  fine  alluvial  materials. 

36.  The  bed  surface  would  be  plane  for  currents  less  than  some 
critical  velocity.  For  currents  slightly  above  critical,  the  first 
features  observed  are  ripples.  Ripples  are  triangular  corrugations  of 
small  height  and  length,  the  upstream  faces  being  flatter  than  the  down- 
stream faces.  At  a higher  current  velocity,  these  irregularities  grow 
in  height  and  in  length  to  form  features  ordinarily  called  dunes.  Both 
ripples  and  dunes  move  downstream.  These  two  features  are  responsible 
for  the  increase  of  resistance  in  an  alluvial  channel.  The  presence  of 
ripples  causes  no  agitation  of  waters  at  the  surface,  while  dunes  gen- 
erate surface  disturbances  in  the  form  of  waves.  This  latter  condition 
has  been  utilized  by  some  investigators  to  formulate  theoretical  rela- 
tions between  the  sand  waves  and  the  surface  waves.  As  the  velocity  is 
increased,  a point  is  reached  where  the  dunes  are  erased  and  the  bed 
becomes  flat.  This  transition  ordinarily  occurs  at  a Froude  number 
less  than  unity  for  ordinary  sands  of  about  1 mm.  As  the  Froude 

1*1 


number  increases  above  unity,  one  may  observe  a sequentially  plane  bed 
with  a smooth  water  surface  and  antidunes.  Antidunes  are  standing  waves 
of  sand  that  are  accompanied  by  corresponding  standing  surface  waves; 
both  the  antidunes  and  the  surface  waves  move  upstream,  grow  in  height, 
and  then  break. 

Initiation  of  sediment  motion 

37.  One  may  proceed  with  dimensional  analysis  to  establish  ex- 
pressions for  the  initial  entrainment  of  sand.  The  quantities  involved 
are:  , the  velocity  of  current  at  the  apices  of  sands;  d , the 

grain  size;  p , the  density  of  water;  pg  , the  density  of  sand; 
v , the  kinematic  viscosity  of  water;  and  g , the  constant  of  gravity. 
Assume  that  the  initial  motion  of  sand  is  an  upward  displacement  due  to 
the  lift  arising  from  the  pressures  around  the  sand  surface  of  the  value 

Fr  = ClPU^d2  (2.1) 

The  constant  in  general  is  a function  of  placement,  of  the  shape  of 

sands,  and  Reynolds  number  U^d/v  . The  uplifting  force  has  to  overcome 
the  weight  of  sand  in  water. 

Wg  = Z (ps  - p)sd3 

Equating  these  two  forces 

pU2  U d 

r r 

± 

(ps  - p )gd  91  v 

concealing  in  the  functional  the  constants  relating  to  placement  and 
shape  of  sands.  In  regular  channels,  the  velocities  may  be  inferred 
from  the  shear  velocity  U#  , U#  = A/p  , where  t is  the  shear  acting 
on  the  channel  floor  bearing  on  the  sands.  Whether  the  flow  be  laminar 
or  turbulent 

u#y 

U = U*  *2  ~ 
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where  U is  the  velocity  at  the  point  of  distance  y from  the  channel 
floor . Thus 

u#d 

U = UM<Ju  

r *2  v 


and  introducing  this  in  Equation  2.3,  the  result  is 


M 
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(p  - p )gd 
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This  initially  was  given  by  Shields.  Shields  has  shown  the  validity  of 

this  expression  using  the  data  of  his  experiments  with  grains  of  varying 

12 

densities  together  with  the  data  from  Casey,  Krammer,  WES,  and  Gilbert. 

The  solid  curve  shown  in  Figure  12  and  referred  to  as  Shields'  curve 

13  l1* 

is  a good  representation  from  Vanoni.  The  curve  was  drawn  by  Rouse 
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Figure  12.  Criteria  of  sand  entrainment  in  terms  of 

shear  stress 
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and  covers  the  range  of  U#d/v  from  2 to  500.  Since  the  boundary  layer 
thickness  6 can  be  expressed  by 

U#6 

— =11-6  (2.5) 

then  the  range  of  d/6  covered  by  the  curve  is  from  0.2  to  50.  The 
minimum  value  of  Shields'  stress  parameter  is 


7F 
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T 

- P )gd 


0.035 


(2.6) 


and  this  at  a point  U#d/v  = 10  and  where  the  sand  diameter  equals  the 
laminar  boundary  layer  thickness.  For  grains  increasing  in  diameter, 
the  parameter  approaches  asymtotically  to  the  maximum  value 


T 

(Ps  - p)gd 


0.06 


(2.7) 


The  extension  of  Shields'  curve  shown  with  a broken  line  in  Figure  12 
was  drawn  by  Vanoni.  It  relates  to  the  condition  where  d/6  is  less 
than  0.2.  This  curve,  having  the  slope  of  -1,  admits  the  equation 


In  Shields'  experiments,  the  critical  stress  is  determined  from  the 
intersection  of  two  lines,  one  corresponding  to  the  condition  of  repose 
and  the  other  to  the  condition  of  increased  particle  motion.  In  the 
data  of  Krammer  and  WES,  critical  stress  is  associated  with  "weak  motion" 
where  few  particles  are  in  motion  which  can  be  counted.  Agreement  of 
data  from  various  investigators  not  shown  in  Figure  12  is  remarkably 
good. 

38.  It  would  appear  that  the  critical  stress  of  sediment  motion 
depends  on  the  state  of  turbulence  in  a cross  section.  Shields  and  the 
other  researchers  studied  channels  where  the  turbulence  of  the  cross 
section  is  fully  developed.  In  contrast  with  this,  White1^  examined 
the  critical  stresses  for  several  sediments  in  the  laminar  and 


kk 


turbulent  boundary  layers  associated  with  the  flow  created  in  a nozzle. 

Contrary  to  Shields'  results,  the  critical  stress  parameter  for  the 

turbulent  boundary  layer  was  independent  of  U#d/v  and  had  the  uniform 

value  of  about  0.10.  This  is  twice  as  large  as  that  read  from  Shields' 

curve,  and  the  enlarged  value  is  attributed  to  the  turbulent  velocity 

fluctuations  which  at  the  bed  may  amount  to  as  much  as  half  the  value  of 

the  local  mean  velocity.  The  stress  parameter  for  the  laminar  boundary 

layer,  U#d/v  less  than  3.5,  was  also  uniform  having  the  valu..  of  0.18. 

13 

This  matter  was  reexamined  by  Vanoni  also  in  a channel  of  developing 
boundary  layer,  using  fine  sediments,  sand,  and  glass  beads.  The  re- 
sults for  the  turbulent  boundary  layer  are  shown  in  Figure  12.  Vanoni 
data  are  in  good  agreement  with  the  extrapolated  Shields'  curve.  This 
may  mean  that  for  U#d/v  less  than  2,  results  regarding  the  critical 
stress  are  independent  of  flow  conditions.  Further  researches  concern- 
ing the  point  are  necessary. 

39-  For  some  problems  such  as  the  stability  of  rocks  of  river 
embankments,  it  would  be  most  helpful  to  turn  to  the  relation.  Equa- 
tion 2.3,  with  the  function  / determined  as  a function  of  U d/v  . 

r 

This  becomes  necessary  in  those  conditions  where  in  principle 
cannot  be  inferred  either  from  the  mean  channel  current  velocity  U or 
from  the  stress  t on  the  rigid  boundaries.  In  the  presence  of  eddies, 
for  example,  takes  on  new  values  and  one  must  now  account  for  the 

increases  realized. 

Valembois  grain  parameter 

1+0.  Valembois  introduces  the  dimensonless  number  G which  is  the 
2 

ratio  of  (U#d/v)  and  t /(p  - p)gd 

0 s 


pv 


This  is  termed  the  grain  parameter,  which  has  a special  significance 

l6 

for  the  transport  of  sediments.  The  Shields'  relation,  putting 
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now  may  be  written  as 


G - <J>^(R#)R# 


or  simply 


I 1 


! ! 


k 


b 

Sr 

§ 


G = 4>2(R*) 


17  ,, 


However,  as  suggested  by  Bonnefille,  it  would  be  more  appropriate 
to  represent  the  grain  parameter  by  a new  symbol  d# 


<1*  = Gl/3  = 


(Pg  - p)g 


pv 


1/3 


(2.9) 


where  d now  enters  as  a first  power.  Hence  the  condition  for  entrain- 
ment could  be  expressed  by  the  relation 


(2.10) 


Bonnefille  noted  that  different  formulae  proposed  by  various  investi- 
gators could  be  reduced  into  the  above  form.  Plotting  results  from 
various  authors,  this  yields  a well-defined  curve  (Figure  13)  with  a 
discontinuity  of  slope  at  R#  = 12  . At  this  value,  the  thickness 
of  the  boundary  layer  and  the  diameter  of  grains  are  of  equal  value 
and  there  is  a separation  between  two  different  hydraulic  regimes. 


For  R#  < 12 


d#  = 2.5Rt/5 


(2.11) 


For  R#  > 12  , the  experimental  points  align  themselves  about  a curve 
with  the  equation 


d#  = 3.8R^8 
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(2.12) 
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Figure  13.  Criteria  of  sand  entrainment  in  terms  of  sand  parameter 
but  it  may  be  said  that  the  law 

r\  /n 

(2.13) 


d#  = 2.7R #/3 


represents  the  phenomena  also  equally  well.  These  expressions  cover  a 
range  of  R#  equal  to  0.2  up  to  10,000,  and  of  d#  from  0.8  up  to 
6,000,  involving  grain  sizes  from  0.04  mm  to  7.5  cm.  From  Equation  2.3, 
introducing  R#  and  d#  , we  obtain 

d#  = 2.56R^3 

which  is  in  agreement  with  Equation  2.11,  and  from  Equation  2.8, 


d*  = 2.08R# 
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(2.14) 
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which  does  not  vary  appreciably  from  Equation  2.13,  the  law  for  R#  of 
the  order  of  unity.  The  utility  of  expressing  the  law  of  sand  entrain- 
ment in  terms  of  d#  , the  grain  parameter,  instead  of  the  Shields' 
stress  parameter  is  that  the  critical  stress  for  a given  d , v , and 
Ps  may  be  obtained  without  trial. 

Critical  current 
velocity  of  sediment  motion 

1+1.  It  appears  that  the  law  of  initial  sand  entrainment  is  well 
established  in  terms  of  the  critical  shear  stress.  On  certain  occasions 
it  may  be  desirable  to  relate  the  entrainment  to  a critical  mean  ve- 
locity of  the  canal.  Resorting  to  dimensional  analysis,  the  general 
expression  for  the  critical  velocity  U is 


(2.15) 


This  functional  relation  may  be  determined  from  entrainment  laws  in- 
volving critical  shear  stresses.  The  relation  varies  depending  on 
the  hydraulic  characteristics  of  the  flow.  The  flow  may  be  hydrodynam- 
ically  smooth  at  the  transition  stage,  or  turbulent  at  the  initiation  of 
sand  movement.  It  must  be  assumed  that  the  Karman  constant  k equals 
0.1+  at  the  start  of  motion.  When  the  grain  diameter  is  much  smaller 
than  the  viscous  sublayer  thickness,  U#d/v  "3.3 
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0.193 


(1.1*1  bis) 


Equation  2.8  may  be  written  as 


Eliminating  the  shear  velocity  U#  between  these  latter  two  equations, 
one  has 


(2.16) 
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which  is  the  critical  velocity  for  the  initial  motion  of  granular  mate- 
rial with  diameters  smaller  than  the  sublayer  thickness.  Note  that  the 

criterion  is  independent  of  the  size  of  the  sediment  particles  when  the 

grain  size  is  small  in  comparison  with  the  viscous  sublayer  thickness. 
Referring  once  more  to  Shields'  curve  in  Figure  12,  it  is  seen  that  the 

stress  parameter  is  minimum  around  R#  = 10  . At  such  regimes  the  sedi- 

ments project  from  the  laminar  sublayer  and  the  shear  velocity  may  be 
approximated  by 

/ u \~l/6 

U*  = 0.123  (j; ) U 

This  is  really  more  valid  for  flow  regime  with  R#  = 60  , and  the  co- 
efficient of  the  formula  needs  to  be  modified  slightly  in  order  that 
the  formula  apply  to  the  transition  from  smoothness  to  complete  rough- 
ness. At  the  moment  we  shall  ignore  the  necessary  corrections.  For 
the  region  in  the  neighborhood  of  R#  = 10  , Shields ' curve  gives 


U#  = 0.030^-  - 1^  gd 


(2. IT) 


Eliminating  U*  between  these  last  two  equations,  one  has 


For  ordinary  sand  p =2.65  and  the  relation  becomes 

s 


1/3 


(2.18) 


Had  we  used  Equation  2.7  applicable  for  gravels  one  would  have  obtained 


(2.19) 


For  ordinary  gravels  pg  = 2.65  and  this  relation  becomes 
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The  Fort  Collins  data  on  the  initial  motions  of  the  sands  are  shown  in 

Table  11.  The  same  table  shows  the  values  of  the  parameters  in  Equation 

2.18.  The  ratios  of  these  parameters  are  plotted  against  U#d/v  in 

Figure  l4  and  indicate  that  the  value  of  the  ratio  of  the  Froude  number 

to  (d/h)1^  tends  to  1.75  in  agreement  with  Equation  2.l8  where  the 

ratio  in  question  has  the  value  1.8l.  This  agreement,  if  not  fortuitous, 

is  rather  surprising.  Equation  2.l8  relates  to  uniform  sand,  whereas 

the  Collins  data  are  for  graded  sand  with  o practically  close  to  1.75. 
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Figure  lU.  Fort  Collins  data  on  critical  velocity 
for  sand  entrainment 


Criterion  of  Formation  of  Bed  Features 


k2.  Since  various  bed  features  affect  the  flow  resistance  in  an 

alluvial  channel  in  a decisive  manner,  it  is  important  that  a criterion 

for  their  formation  be  derived.  The  problem  was  approached  first  by 
19 

Liu.  Noting  that  at  the  instant  of  ripple  formation,  the  water  sur- 
face in  the  channel  is  smooth  and  without  waves,  an  instability  of  flow 
near  the  bottom  was  considered  to  be  the  cause  of  the  ripples.  The 
analysis  proceeds  analogously  as  with  the  instability  of  an  interface 
between  two  liquids  of  different  density.  The  conclusion  is  that 


where  W is  the  fall  velocity.  Now,  obviously. 


and  hence  expressing  this  in  a dimensionless  form 


w2  * [(t)-  1 *d]  *M) 


Introducing  this  in  Equation  2.21 


_ (u*d) 

!p  „/p  - l)gd  ^ ' v / 


(2.22) 


The  conclusion  is  that  if  Equation  2.21  is  valid  to  represent  the  initial 
formation  of  ripples  Equation  2.22  is  also  valid.  Liu  employed  both 
equations  for  the  examination  of  data  of  his  experiments  and  those  of 
others  and  remarks  that  there  is  better  correlation  with  the  parameter 
U#/W.  The  scatter  of  data  is  due  probably  to  the  uncertain  values  of 
U#  , which  would  be  twice  as  large  in  Equation  2.22.  We  shall  pursue 


our  examination  of  the  criterion  of  formation  of  bed  features,  applying 

Q 

Equation  2.22  to  the  Fort  Collins  data. 


U3.  The  method  followed  to  determine  the  shear  velocities  Ik 


associated  with  the  initiation  of  the  various  bed  feature  formations  is 
as  follows.  For  a specified  sand,  the  shear  velocities  U#  and  the 
corresponding  sediment  transport  quantity  q^.  are  read  from  the  data 
tables  given  in  Reference  18  and  are  plotted  one  against  the  other  on 
semilogaritbmic  coordinate  paper,  as  shown  in  Figure  15.  The  points  be- 
longing to  different  bed  features  are  differently  indicated.  Fortunately, 
the  laws  giving  the  dependence  of  q^.  on  U#  are  different  for  the 
different  bed  features,  and  this  allows  us  to  single  out  the  start  of 
the  various  bed  features.  The  open  squares  show  the  values  of  U#  cor- 
responding to  the  initial  motion  of  sediments.  The  intersection  of  the 


straight  line  passing  through  the  solid  squares,  the  marks  for  ripples, 

with  the  (i  = 1 x 10  ^ axis  is  taken  at  the  start  of  ripples.  Thus, 

~ —5  2 

for  sand  d = 1.1+7  * 10  ft  and  for  viscosity  v = 1.25  * 10  ft  /sec. 


U#  = 0.075  ft/sec  is  the  shear  velocity  at  the  start  of  ripple  forma- 


tion, U„  = 0.112  ft/sec  at  the  start  of  dune  formation,  and  U#  = 0.22 
ft/sec  for  transition.  This  and  the  corresponding  data  for  the  remaining 
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Figure  15.  Shear  velocity  versus  sand  transport  in  Fort  Collins 

8-ft  flume 

sands  are  given  in  Table  12.  In  the  Fort  Collins  tests,  the  shear 
stress  is  determined  as  t = pgHS  . This  means  that  sidewall  frictions 
are  neglected  and  corrections  are  needed  to  obtain  the  stress  of  the 
bottom. 

Sidewall  corrections 

44.  The  Fort  Collins  tests  were  carried  out  in  two  different 
flumes.  In  one,  with  the  channel  width  of  8 ft,  the  maximum  value  of 
the  ratio  of  water  depth  to  channel  width  was  about  0.25.  The  correc- 
tions in  the  determination  of  the  channel  bottom  shear  stresses  would  be 
small.  In  the  other,  with  the  channel  width  of  2 ft,  the  maximum  value 
of  the  ratio  of  water  depth  to  channel  width  was  about  0.50.  For  the 
latter,  sidewall  corrections  could  be  considerable.  Ordinarily,  two 
methods  are  used  to  correct  for  wall  friction:  one  is  the  Johnson 
method,  and  the  derivations  are  based  on  the  assumption  of  a logarithmic 
velocity  distribution  of  smooth  walls;  the  other  is  the  Einstein  method 
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Table  12 


Data  on  the  Start  of  Ripples,  Dunes,  and  Transition 
Deduced  from  Fort  Collins  Studies 


Sand 

d50  x io3 
ft 

Start  of 

Ripples 

Start 

of  Dunes 

Transition 

u* 

v x io5 

u* 

v x io5 

u# 

v x io5 

0.622 

0.065 

1.15 

0.11+ 

1.15 

0.185 

1.15 

0.880 

0.060 

1.20 

0.15 

1.13 

0.19 

1.15 

0.910 

0.060 

1.38 

0.15 

1.15 

0.19 

1.21 

1.U7 

0.075 

1.28 

0.12 

1.25 

0.22 

1.25 

1.53 

0.090 

1.20 

0.13 

1.20 

0.26 

1.20 

3.05 

— 

— 

0.11 

1.09 

0.33 

1.09 

1.05 

0.062 

1.10 

0.135 

1.10 

0.20 

1.10 

1.08* 

0.080 

1.08 

0.129 

1.08 

0.23 

1.08 

1.08** 

0.068 

1.12 

0.115 

1.03 

0.20 

1.03 

l.?6 

0.065 

1.20 

0.100 

1.10 

0.31 

1.08 

* Uniform. 

**  Graded . 

and  assumes  a Manning's  form  of  velocity  distribution.  Here  the  method 
of  Einstein  will  be  used,  but  first  it  is  desirable  to  establish  the  de- 
sired formula  using  the  Johnson  approach.  The  channel  cross  section  Aq 
may  be  divided  into  two  parts,  and  A^  : the  first  is  the  domain 

of  influence  of  the  frictional  action  of  the  two  walls  and  the  second, 
the  domain  of  influence  of  the  frictional  action  of  the  bottom.  Assume 
that  at  the  points  of  the  lines  separating  the  two  areas  the  stresses 
vanish.  This  assumption  is  correct  as  an  approximation.  To  be  abso- 
lutely valid,  the  velocity  rates  normal  to  the  lines  of  separation  should 
vanish,  which  is  not  always  the  case.  The  second  approximation  is  that 
the  average  values  of  the  velocities  in  the  two  regions  A^  and  A^  are 
equal  to  each  other.  This  is  also  only  approximately  correct,  as  the 
mean  velocity  in  A^  is  smaller  than  that  in  A^  . Subject  to  this 


5*+ 
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resolution  of  forces  in  a channel  section  gives  the  relations 


assumption 


S = slope  of  energy  gradient 
p = density  of  vater 

tq  = average  shear  stress  of  all  the  rigid  boundaries 
= average  shear  stress  of  the  bottom 
= average  shear  stress  of  the  two  walls 
The  channel  wetted  perimeter  is  p^  ; p^  and  pw  are  the  wetted  lines 
of  the  sides  and  of  the  bottom,  respectively.  Let  B be  the  channel 
width  and  H the  depth;  thus 


then 


Introducing  the  definitions  of  the  coefficient  of  resistance  X 


where  U is  the  mean  velocity  of  the  overall  cross  section  A or  the 

o ’ 

subsections  A^  and  A , Equation  2.23  becomes 


X p 
w *w 


A p 
o 


(2.24) 


Ik  _ i + 2H  2H 

To  B " B A 

° o 


Thus,  the  ratio  of  the  effective  shear  stress  on  the  bottom  to  the  stress 

of  the  entire  channel  may  be  obtained  if  the  ratio  A /A  is  known.  The 

o w 

shear  stresses  in  terms  of  the  hydraulic  radii  are 


To  = pRogS 


Tb  = PVS 


Tw  = PVS 


Comparing  these  with  the  expressions  of  the  stresses  in  terms  of  the 
current,  given  above,  we  have 


R A R A 

_w  - b b 

Rn  " A and  r"  = — 
o o o o 


(2.25) 


The  coefficient  of  resistance  A in  terms  of  Manning's  n is 


2.208R 

where  R is  measured  in  feet  and  g is  32.2  feet  per  second  per 
second.  Hence 


A = 29.1n2R-1/3 


1 


specializing 


X = 29. ln2R“1/3 
w w w 


Substituting  from  Equation  2.25  for  R^  and  solving  for  X^ 


X = 12 

v 


•>*"&) 


and  then 


~ = 12.5n3/2X"3/4R"1/1+ 
X woo 

o 


(2.26) 


In  channels  where  the  sidewalls  are  constructed  from  glass  or  aluminum 
an  appropriate  value  n^  is  0.01.  For  such  channels,  using  Equation 
2.26  in  connection  with  Equation  2.24 


r ■ 1 * t?  - (o.oissx^V1^)  f 


(2.27) 


Introducing  the  Darcy-Weisbach  coefficient  f , fQ  = 4Xq 


r ■ 1 (°-0353r3/VlA)  & 


(2.28) 


In  the  data  of  the  Fort  Collins  tests  the  stresses  t and  p are 


computed  as 


no  , „ 8rHS 

t = pgHS  and  f = SL„- 

U 


(2h 
i + — ) 
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Introducing  these  in  Equation  2.28 


= (l  - 0.0353f-3/Sr1/U 


(2.29) 


The  values  of  H , B , f , and  r are  read  from  the  data  tables  of 

^ Q 

the  Fort  Collins  tests.  More  generally  it  may  be  desirable  to  resort 
to  Manning's  nQ  of  the  composite  channel,  n^  of  the  walls,  and 

of  the  bottom.  From  the  analysis  and  omitting  the  details  of  the 
derivations 


R /n  \3,/2  R /rn  \3/ 

R \n  / ’ R = \n  I 

o \ o/  o \ o/ 

3 

— = (— \ = i + 2g.  2H  /M 

t l R ) B ‘ Bln 

o \ o / \ o / 


(2.2^  bis) 


Usually  it  is  the  custom  when  analyzing  the  characteristics  of  bed 
features  observed  in  the  laboratory  composite  channel,  to  replace  the 
channel  with  an  infinitely  wide  channel  of  energy  gradient  Sq  , cur- 
rent velocity  U , and  depth  . It  is  equally  serviceable  to  choose 

the  wide  channel  with  the  energy  gradient  S,  (=S  R /R,  ) , current  veloc- 

b o o b 

ity  U , and  the  depth  H . In  the  studies  of  the  bed  features  sub- 
sequently adopted  the  latter  selection  will  be  followed. 

First  appearances  of 
ripples  and  dunes  according 
to  the  Fort  Collins  data 

^5.  The  formula  of  correction  Equation  2.29  was  applied  to  the 

data  shown  in  Table  12.  In  general,  during  the  regimes  of  ripples  or 

dimes,  according  to  these  computations  bottom  shear  stress  r,  is 

b 

smaller  than  t , the  channel  shear  stress  evaluated  by  neglecting  the 
resistance  of  walls  by  7 percent  in  the  8-ft-wide  channel.  The  cor- 
responding reduction  in  the  2-ft-wide  channel  is  about  25  percent. 

Data  of  Table  13  incorporate  these  corrections  and  then  represent  the 
shear  velocity  proper  of  the  channel  bed;  the  plotting  in  Figure  l6  is 
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Table  13 


Criteria  of  Appearance  of  Ripples,  Dunes,  and 
Transition  with  Sidewall  Corrections 


F#  = U#/(p  /p  - l)gd 


R»  = U.d/v 


Sand 

Start  of  Ri] 

spies 

Start  of 

Dunes 

Transition 

d50  - 103 

V2 

R* 

p2 

r# 

R# 

F2 

r * 

R# 

0.62 

0.120 

3.1+ 

0.552 

5.9 

0. 842 

9.63 

0.88 

0.073 

4.2 

0.460 

11.3 

0.705 

14.0 

0.91 

0.070 

4.4 

0.435 

11.5 

0.670 

13.8 

1.1+7 

0.067 

8.6 

0.173 

13.5 

0.566 

25.0 

1.53 

0.093 

11.0 

0.192 

16.0 

0.760 

32.2 

3.05 

— 

— 

0.065 

29.8 

0.612 

89.0 

1.05 

0.059 

5.4 

0.268 

11.4 

0.590 

17.4 

1.08* 

0.096 

7.4 

0.234 

11.5 

0.782 

21.6 

1.08** 

0.051 

6.0 

0.186 

10.7 

0.590 

18.8 

1.76 

0.038 

8.9 

0.097 

10.6 

0.880 

46.9 

* Uniform. 

**  Graded. 

from  this  table.  Despite  the  scatter  of  the  plotted  points,  straight 
lines  should  suffice  to  represent  them  as  a first  approximation.  The 
line  for  the  ripple  formation  cuts  the  Shields'  curve  at  a point 


7 rr-r  = 0.032,  = 12.0 

(s  - l)gd  v 

and  accordingly 


(ps/p  - l)gd  ^ (12.p)2  _ 


0.032 


= 1+500 


d„=  '■'Jo : 


l)g 


d = 16.5 


(2.30) 


f 


- ^ » , < 


Figure  1 6.  Domains  of  the  bad  features  in  the  Fort  Collins  tests 

which  is  the  criterion  of  the  limiting  grain  size  for  the  ripple  forma- 
tion. For  quartz  sand,  this  limiting  value  is  0.65  mm  on  the  basis 
that  the  viscosity  of  water  is  v = 1.0  * 10  ft  /sec.  With  the 
sands  having  grain  sizes  larger  than  this,  no  ripple  would  be  formed. 
For  coal  particles,  Ps/p  = 1-30  , the  corresponding  limiting  value 
is  1.14  mm.  Reductions  shown  in  Figure  16  are  mostly  from  Fort 
Collins  data  for  sands  with  sizes  smaller  than  0. 6k  mm,  excepting  one 
dj-Q  = 0.93  mm.  In  the  tests  with  this  latter  sand  no  ripples  were 
formed;  the  surfaces  remained  smooth  with  sand  transportation  present. 
Thus  information  regarding  the  formation  of  bed  features  that  might 
be  abstracted  from  Figure  l6  need  be  confined  to  sands  with  grain 
sizes  smaller  than  the  criterion  d#  < 16.5  • The  straight  line  drawn 
through  the  points  that  mark  the  beginning  of  ripples  has  the  slope 
of  -1.  Accordingly,  the  channel  bottom  shear  velocity  at  the  start  of 
the  ripple  formation  is  expressed  by  equation  of  the  line  of  the  start 
of  ripples 
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(2.31) 


U#  = 0.72[(s  - l)gv]1/3 

The  slope  of  the  straight  line  drawn  through  the  points  that  mark  the 
beginning  of  the  formation  of  dunes  is  equal  to  -3.  Therefore,  channel 
bottom  shear  velocity  at  the  start  of  the  dunes  is 

U*  = 3.06[(s  - l)gd]l/5(!)  (2.32) 

The  expression  of  the  criteria  for  the  start  of  the  ripple  formation  is 
independent  of  the  size  of  the  grains.  In  contrast  with  this,  the 
criteria  for  the  start  of  the  dune  formation  does  depend  on  the  size  of 
grains.  Because  of  the  scatter  of  points,  the  true  form  of  the  cri- 
teria by  the  present  method  cannot  be  obtained.  The  transition  from 
dunes  to  flat  bed  occurs  at  the  shear  velocity. 

U#  = 0.86V(s  - l)gd  (2.33) 

Note  that  the  shear  stress  is  independent  of  viscosity.  The  intersec- 
tion of  the  curves  marking  the  start  of  the  formation  of  dunes  and 
transition  yields  the  critical  of  value  of  d#  for  the  condition  that 
ripples  are  followed  by  transitions.  To  obtain  the  intersection,  we 
consider  the  short  dashed  line  passing  through  the  point  corresponding 

to  d =0.19  mm.  This  line  has  the  same  slope  as  the  full  line  giving 
m 

the  start  of  the  dune  formation  for  most  of  the  observed  data.  The 
intersection  yields 


or 


U*d 


T__7_r  , o.75;  — = 1.9 


(s  - l)gd3 

2 

v 


2b.  5 

0.75 


32.8 


(2.3b) 


cL#  = 3.2 


6l 


f 


I! 
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This  is  the  critical  value  for  the  condition  that  ripple  is  followed  by- 
transition  without  the  intermediate  dune  formation. 

b6.  The  field  of  validity  of  the  above  relations  is  the  region 
between  the  lines  d#  = 3.2  and  d#  = l6.5  in  Figure  16  and  this  is 
the  field  proper  for  sand  of  size  less  than  the  critical,  with  its 
dimensionless  value  as  d#  = l6. 5 . The  line  d#  = 5.0  goes  with  the 
sand  of  d^^  = 0.19  mm  of  the  Fort  Collins  studies,  the  smallest  sand 
used  in  these  studies.  The  marked  data  deviations  noted  are  not  en- 
tirely due  to  observational  errors  but  are  brought  about  by  the  granulo- 
metric dissimilar  size  distribution.  For  example,  the  sand  in  the 
series  referred  to  as  of  median  size  of  1*7  mm  was  produced  from  l*5-mm 
sand.  Coarser  sand  was  retained  and  the  fines  were  washed  away  in 
overflow.  Although  as  a result  the  change  in  gradation  was  small,  from 
a = 1.60  to  a = 0.5**  , with  the  fines  removed  the  starting  of  the 
ripples  was  markedly  delayed.  Precise  determination  of  the  criteria  of 
the  appearances  of  bed  features  cannot  be  made,  apparently,  if  the 
effects  of  gradation  are  ignored. 

The  Chatou  studies  on  the 
appearance  of  bed  features 

1*7.  In  the  above,  the  question  of  the  appearance  of  bed  features 
for  sands  with  sizes  smaller  than  the  critical,  d#  < 16.5  , was  con- 
sidered. Projection  from  this  to  the  cases  where  the  sand  grain  size 


is  greater  than  the  critical  is  not  allowed.  Fortunately,  for  the 

desired  information  one  may  now  turn  to  the  Chatou  studies  carried  out 

20 

by  Chabert  and  Chauvin.  The  experiments  were  made  in  an  inclined 
flume  with  glass  walls  21  m long,  80  cm  wide,  and  70  cm  deep  with  a 
movable  bed  of  15-cm  thickness.  Eleven  different  materials  were  con- 
sidered with  densities  ranging  from  1.075  to  2.62  and  with  particles 
with  sizes  ranging  from  0.17  mm  to  3.0  mm.  Seven  of  the  materials  had 
sizes  greater  than  0.55  mm  and  with  these,  ripple  formation  was  not 
observed.  The  authors'  description  of  the  behavior  of  these  materials 
under  the  action  of  increasing  tractive  force  is  paraphrased  as  follows 
Following  the  initial  entrainment  of  the  sediment,  material  is  trans- 
ported without  undulations.  The  flats  are  the  only  possible  forms  of 


the  bed.  If  a small  swelling  is  artificially  produced,  it  is  soon 
wiped  out.  Very  long  flat  plateau  waves  are  noted  which  eventually 
evolve  into  dunes  of  large  height,  closely  spaced,  and  often  covered 
with  smaller  dunes.  Omitting  the  experimental  points  of  the  observa- 
tions, the  appearance  of  bed  forms  is  very  neatly  represented  by  the 
curves  in  Figure  17.  Delineations  are  indicated  separating  the  various 


Figure  17.  Domains  of  ripples  and  dunes  as  shown  by 
Chabert  and  Chauvin 


regimes  from  each  other.  The  curve  for  the  entrainment  of  the  particles 
is  somewhat  displaced  downward  from  the  curve  of  Shields.  The  area  to 
the  left  of  the  line  d#  = 15  represents  the  domain  of  the  ripples. 
According  to  that,  the  particles  with  the  dimensionless  diameter 


d*  > 15  will  not  engender  ripples  but  flats.  This  is  shown  in  the 
figure.  The  curve  that  separates  the  domain  of  dunes  from  the  domain 
of  flat  bottoms  to  the  right  of  the  line  d*  = 15  shows  that  dures 
will  appear  at  the  instant  that  the  bottom  shear  velocity  reaches  the 
average  value  of: 


(pg/p  - l)g<i 


= 0.055 


(2.35) 


The  data  of  Figure  l8  are  in  agreement  with  some  of  the  findings  shown 
here.  Previously,  from  Figure  l6,  it  was  deduced  that  the  critical 
diameter  of  sediment  particles  associated  with  ripple  formation  should 


Figure  l8.  Friction  loss  due  to  channel  bed  features 
in  alluvial  streams 

be  d*  = 16.5  • This  compares  well  with  the  criterion  d*  = 15*2  of 
figure  established  more  definitely.  In  the  Fort  Collins  study,  no 
ripples  were  formed  with  sand  of  0.93  mm  and  of  dimensionless  size 
d*  = 23.2  . According  to  Table  13  and  Figure  l6,  when  dunes  first 
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appeared  the  shear  velocity  for  this  sand  was 


(ps/p  - l)gd 


= 0.065 


a value  which  is  in  close  agreement  with  the  one  shown  above. 


I 


If 


Resistance  of  Alluvial  Channels 


Einstein-Barbarossa  analysis 
of  alluvial  channel  resistance 

48.  It  is  postulated  in  this  analysis  that  the  total  resistance 

t of  an  alluvial  channel  may  be  broken  into  two  parts  x^  and  x^  , 

the  first  being  due  to  friction  of  granular  surface  and  the  second  to 
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the  form  resistance  of  the  ripples  and  dunes.  Therefore 


x X1  + t2 


(2.36) 


Dissipation  of  energy  associated  with  these  two  frictional  stresses 
occurs  in  two  areas,  and  . These  areas  are  not  separate  but 

are  interlocked  and  have  the  same  wetted  perimeter  p as  the  entire 
cross  section  A . Correspondingly,  will  be  the  hydraulic  radius 

associated  with  the  granular  resistance  and  R^  that  associated  with 
the  form  resistance;  that  is 


Rx  = A1/P 


R2  = Ag/p 


Denoting  the  hydraulic  radius  of  the  entire  channel  R and  as 
A±  + A2  = A , 

h = R-^  + R2 


(2.37) 


The  current  mean  velocity  in  relation  to  the  granular  surface  resistance 
is 
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(2.38) 


i 

j 


- 5-75  log  ^12,21  -tj 

U#1  =VTl/p  =VgR  Is  (2.39) 

where 

S = slope  of  the  energy  line 

A = an  effective  roughness  height  related  to  the  height  of  the 
actual  asperities  k 

Where  the  bed  material  is  not  uniform  k is  taken  to  equal  dg^  , the 
grain  size  65  percent  finer  by  weight.  Einstein's  sediment  transport 
theory  assumes  that  only  the  stress  action  on  the  grains  is  responsible 
for  the  movement  of  bed  material  and  the  subsequent  formation  of  bed 
features.  The  stress  parameter  pertinent  in  the  transportation  is  the 
ratio  of  submerged  weight  of  a particle  to  the  shear  stress  action  on  it. 


(Ps  - p)gd3 
pgR  Sd2 


(s  - l)d 
R1S 


(2.1+0) 


For  quartz  sand  s = 2.65  and  if  the  bed  material  is  not  uniform  in 
size,  d is  taken  to  equal  d^  . With 

U*2  =VT2/,p  =V®R  2S  (2.1+1) 


the  form  resistance  of  bed  is  assumed  to  depend  on  ip'  and  thus 


* 


(2.1+2) 


The  form  of  the  latter  function  was  determined  by  Einstein  and 
Barbarossa  on  the  basis  of  observation  from  a few  American  alluvial 
rivers  with  grain  sizes  ranging  from  0.19  mm  to  U. 30  mm.  Average  values 
of  these  determinations  are  shown  by  circles  in  Figure  18.  The  curve 
in  dashes  has  been  referred  to  as  the  Einstein-Barbarossa  bar  curve. 

We  have  drawn  the  straight  line  in  full  in  the  figure  as  an  approximate 


representation  to  the  observed  values.  This  has  the  equation 

,-1/2 


U#2  p R1S/ 


(2.43) 


which  will  be  used  subsequently. 

49.  Expressing  the  quantity  A appearing  in  Equation  2.38  as 


A = mk 


(2.44) 


the  proportionality  factor  will  depend  on  the  circumstances  whether  the 
channel  is  hydrodynamically  smooth  or  rough  in  transition.  When 
U#1k/v  is  less  than  3.5,  the  channel  is  hydrodynamically  smooth  and 


™ - 3.3v  . _ 3.3v 

m = - — p or  A = — 

u«!k  U#1 


(2.45) 


when  U^k/v  is  greater  than  60  or  70,  the  channel  is  rough  and 


m = 1 or  A = k 


U k 

3.5  < — — < 60 


(2.46) 


m will  be  a numeric  as  shown  in  Table  l4.  The  determination  of 
may  be  carried  out  in  the  following  manner.  Equation  2.38  may  be 
written  as 


■^2-.  5.75  (^) 


(2.47) 


which  determines  U//gAS  as  function  of  R^/A  . In  Table  15  the  tabu- 
lar relation  between  these  two  quantities  is  given.  The  determination 
of  R^  when  U , R , k , and  S are  given  is  made  by  trial.  Assume 
a value  of  R^  somewhat  smaller  than  R . Evaluate  U#1  using  Equa- 
tion 2.39*  Establish  n to  obtain  A . Compute  U//gAS  and  find  the 
corresponding  value  R^/A  using  the  entries  of  Table  15.  From  this 
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channel  roughness  is  investigated.  The  materials  used  were  fine  sand 
of  median  diameter  0.23  mm,  coarse  sand  of  median  diameter  4.6  mm,  and 
coal  of  median  diameter  2.2  mm.  The  data  of  resistance  of  the  tests, 
with  sidewall  effects  removed,  were  necessitated  to  evaluate  U/U,0  in 
the  manner  suggested  by  the  Einstein-Barbarossa  analysis;  results  are 
shown  in  Figure  18.  The  results  from  the  fine  sand  are  in  good  agreement 
with  the  curve  for  the  alluvial  rivers,  whereas  the  results  for  the 
coarse  sand  and  for  the  coal  fall  well  above  the  curve,  indicating  less 


resistance.  With  the  fine  sand  tests,  ripples  were  present;  with  the 
tests  using  coarser  particles,  sand  or  coal,  the  channel  bottom  was  flat. 
To  examine  if  these  conditions  are  in  accordance  with  the  criteria  de- 
veloped previously,  the  dimensionless  grain  sizes  d#  are  calculated  in 
Table  l6  for  the  tests  with  coarse  sand  and  coal.  For  the  coarse  sand. 

Table  l6 


Dimensionless  Diameters  of  Franco  Test  Sands 


F2 

r* 

II 

c 

*^ro 

w 

- l)gd;  R* 

— TT  A / w • ^ — K 

s - l)gl 

1/3 

d 

u*a/v,  d#  - 

2 J 

V 

Run 

2 

v ft  /sec 

F2 

R* 

4* 

Bed  Feature 

Coarse  Sand 

: d50=15 

i.4  X 10'3  ft,  s 

-1  = 1. 

68 

36 

0.93 

X icf| 

0.0381 

0.310 

99.0 

Flat  bed 

39 

0.93 

* 10Is 

0.0493 

0.333 

99.0 

and 

4l 

0.93 

x 10J? 

0.0474 

0.327 

99.0 

transport 

37 

1.67 

x 10J? 

0.01+14 

0.170 

88.0 

of 

38 

1.67 

x 10  J? 

0.0488 

0.185 

88.0 

sediment 

1*2 

1.67 

x 10  5 

0.0474 

0.202 

88.0 

Coal : d 

50  * T-8  - 

10'3  ft,  s - 1 = 

° 0.30 

23 

0.93 

x 10J? 

0.0600 

56.3 

37.6 

Flat  bed 

25 

0.93 

x ioj| 

0.0751 

65.0 

37.6 

and 

27 

0.93 

x 10J? 

0.0393 

47.4 

37.6 

transport 

29 

0.93 

x kT: 

0.0870 

70.8 

37.6 

of 

22 

1.67 

x 10"^ 

0.0512 

29.0 

26.0 

sediment 

pli 

1.67 

x 10j 

0.0546 

29.9 

26.0 

26 

1.67 

x 10" ' 

0.0421 

26.2 

26.0 

28 

1.67 

x 10' 5 

0.0782 

35.8 

26.0 

the  dimensionless  grain  sizes  are  99.0  and  88.0  for  two  kinematic  vis- 
cosities; for  the  coal,  these  sizes  are  37-6  and  26.0.  These  values  are 
all  larger  than  the  critical  d#  = 15-6  , and  therefore  ripples  should 
never  be  formed.  Examining  the  table  it  is  seen  that  the  maximum  stress 
parameter  for  the  two  kinematic  viscosities  noted  in  the  entries  for 

p 

coal  is  about  F#  = 0.087  ; that  for  the  two  kinematic  viscosities  for 

2 

coarse  sand  is  about  F#  = 0.01+9  • One  expects  that  the  appearance  of 
the  dunes  for  these  larger  particles  will  be  for  a higher  stress  param- 
eter than  this.  The  Chatou  data,  as  it  will  be  remembered,  places  the 
2 

criteria  at  F#  = 0.055  • The  determination  from  the  Fort  Collins  data 
2 

gave  F*  = 0.065  for  d^Q  = 0.93  mm.  Tentatively,  one  may  place  the 

criterion  for  the  appearance  of  dunes  with  larger  particles,  d#  > 15  » 

P 

as  F#  = 0.69  , the  average  of  the  above  three;  that  is  0.087,  0.055, 
and  0.065. 

51.  The  comparison  of  laboratory  data  with  field  data  in  rela- 

23 

tion  to  \J>'  and  U/U#2  was  previously  made  by  Vanoni  and  Brooks 
using  their  own  laboratory  tests  and  also  those  of  Barton  and  Liu.  In 
the  latter  case,  there  was  good  agreement  with  the  Einstein-Barbarossa 
bar  resistance  curve.  In  their  own  tests,  however,  although  the  agree- 
ment was  satisfactory  for  large  values  of  \p'  , considerable  deviation 
was  observed  with  points  of  small  \|>'  , suggesting  that  with  increasing 

velocity  the  "flume  rivers"  get  smooth  faster  than  natural  rivers.  The 
24  / 

studies  of  Smith  reveal  that  relative  depth  R^/A  is  another  param- 
eter to  consider  in  the  relation 


Although  Smith's  modifications  produce  improvements,  in  our  subsequent 
application  the  simplicity  of  the  Einstein-Barbarossa  procedure  will  be 
preferred. 

Equivalent  sand 
roughness  in  WES  tests 

52.  As  mentioned  previously,  in  the  WES  tests  with  coarse  sand 
and  coal  particles  ripple  formation  was  absent.  Subsequent  to  entrain- 
ment, the  surfaces  of  the  beds  were  smooth.  With  coal  particles,  for 
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1 1 


example,  the  bed  was  flat  with  waves  across  the  width  of  channel  8 to 
10  ft  apart.  Although  there  was  movement  of  bed  material,  the  resis- 
tance to  flow  was  very  small.  Since  there  was  movement  of  bed  material, 
the  question  would  arise  if  this  fact  would  have  an  effect  on  the 


effective  or  equivalent  sand  roughness  of  the  bed  material.  After  a 

25 

careful  study  Vanoni  and  Nomicos 


indicate  that  Harman' s constant 
would  decrease  if  the  material  is  suspended.  As  particle  suspensions 
were  avoided,  in  the  WES  tests  it  is  assumed  that  for  these  cases  the 
Karman  constant  has  the  usual  value  and  the  sand  equivalent  roughness 
may  be  computed  by 


U, 


U, 


5.75  log 


12. 8R 


*1  S 

provided  that  the  particle  sizes  are  larger  than  the  laminar  boundary 
layer  thickness.  For  the  coarser  sand,  the  ratio  k/6  varies  from  15 


to  30  and  with  the  coal  particles,  from  3 to  6.  Computing  kg  on  this 


basis  the  values  of  ks/k  are  shown  in  Figure  19 . The  size  k and 


I 


j 


J 


p 

i 

i 

ft 


i 


the  gradation  number  a of  the  particles  are  also  shown  in  the  figure. 
The  gradation  number  is  computed  by 


The  ratio  k /k  is  plotted  against  q*/UH  , where  q*  is  volume  rate 

5 S S 

of  flow  of  particles  expressed  as  cubic  feet  per  hour  per  foot  of  chan- 
nel width.  Here,  U is  the  mean  current  velocity  in  feet  per  second 
and  H , depth  in  feet.  Although  there  is  some  scatter  in  the  points, 
perhaps  the  line  as  drawn  is  acceptable.  On  this  basis,  one  may  ven- 
ture to  state  that  the  equivalent  sand  roughness  of  particles  when  in 
motion  is  increased.  This  matter  was  further  examined  in  the  corre- 
sponding data  of  Fort  Collins  studies  for  plane  surfaces,  including  also 

the  cases  of  initial  entrainment.  For  these  later  ones,  q*  is  nil. 

s 

Results  of  the  determinations  are  shown  in  Figure  20.  Although,  in 

general  there  is  the  indication  that  k /k  varies  with  the  sediment 

s 

discharge,  the  dependence  is  uncertain.  It  is  also  surprising  to  find 

such  large  values  for  the  ratio  as  10  to  30. 

53.  It  would  thus  appear  that  the  law  of  surface  friction,  in 

the  presence  of  particle  motion,  will  be  different  from  the  law  given 

by  Equation  2.38.  To  account  for  the  difference,  multiply  the  actual 

26 

asperity  size  by  a constant  number.  Engelund,  for  example,  guided  by 
his  own  experiments  suggests  that  d^_  be  multiplied  by  2.  Lovera 
and  Kennedy  studied  this  problem  with  considerable  detail.  Relying 
on  information  mostly  from  natural  alluvial  studies  and  also  partly 
from  laboratory  findings,  they  suggest  that  the  surface  friction 
coefficient 

f,  U*1 

Ai  = ~x = 2 — I (2*W) 

u 

be  expressed  as  a function  of  R/d<-Q  and  UR /v  . From  one  of  their 
figures  where  4X  is  graphed  logarithmically  against  UR/V  for  constant 
values  of  R/d^g  » one  may  derive  the  relation 
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_2 

X = 0.21  ^ (2.1*9) 

50 

2 2 

for  R/d<jo  less  than  50  * 10  but  greater  than  5 x 10  . The  attempt 
to  correlate  the  previously  considered  WES  and  Fort  Collins  data  with 
the  above  equation  was  not  successful.  The  lack  of  conformity  may  be 
noted  also  from  Figure  21  where  X^  is  plotted  against  UR/v  . Each 
symbol  refers  to  a constant  R/d^g  ‘ '^lie  curve  giving  the  lower  bound 
of  the  data  is  the  friction  curve  of  hydrodynamically  smooth  surfaces, 
computed  from  the  Prandtl  formula 

U U*R 

jr-  = 3.25  + 5. 75  log  — (2.38) 

v 

with  individual  values  given  in  Table  IT.  In  the  Lovera  and  Kennedy 
studies,  the  majority  of  points  fall  in  a narrow  band,  with  the  median 
line  of  the  band  of  the  relation 


rr  U#R 

7j  = -3.1*  + 5. 75  log  — 
u*q 

or 

U R 

~~  = -3.*+  - 5-75  log  r~-  + 5.75  log  — (2.50) 

u#q  U#1  v 

which  will  be  referred  to  as  the  modified  Prandtl  formula.  The  tabular 
values  are  given  in  Table  l8.  The  second  curve  in  Figure  17  is  a plot 


of  this  last  equation.  Both  Equations  2.1*9  and  2.50  may  be  used  to 


obtain  X . It  is  not  clear  if  such  evaluations  would  be  valid  for 
-1- 

alluvial  channels  where  R/d^-0  Is  of  the  order  of  1 * 10  or  greater. 
In  these  cases  it  may  be  sufficient  to  determine  X^  on  the  basis  of 
Equation  2.38. 


Coeffi 

cient  of  S 

Table  17 

urface  Friction  for 

Hydrodynami c ally 

Smooth  Surfaces 

^ x lo- 


ll.2 
7.78 
6.16 
5.00 
It. lit 
3.1*6 
2.91+ 


1.05  x ioJ 

2.40 

6.57  k 

1.63  x 10 

4.00 

9-77  , 

2.30  x 10-5 


U/U*! 

A1  x io3 

UR/v 

28 

2.54 

5.22  x 10^ 

30 

2.20 

1.28  x 10° 

32 

1.96 

3.72 

34 

1.63 

7.50 

36 

1.55 

1.73  x 10‘ 

38 

1.38 

4.23 

40 

1.25 

9-95 

Table  18 

Coefficient  of  Surface  Friction  for  Plane  Surfaces 


of  Alluvial  Channels 


X x io- 


13.8 

10.2 

8.88 

7.80 

6.90 

6.l6 

5-54 


5.76  x ioJ 

1.48  x io 

2.24 

3.81 

5.89 

9.34 

1.55  x 10" 


u/un 

X±  x 103 

UR/v 

20 

5.00 

2.40  x 105 

21 

4.52 

3.63 

22 

4.12 

5. 63  6 

24 

3.46 

1.4l  x 10 

26 

2.90 

3.39 

28 

2.54 

8.13  7 

30 

2.20 

1.82  x 10' 

Alam  and  Kennedy  analysis 
of  alluvial  channel  resistance 

54.  In  this  analysis^  the  surface  energy  gradient  SQ(=Sb)  is 
divided  into  two  parts  S1  and  S2  instead  of  dividing  the  hydraulic 

radius  R into  two  parts.  This  is  in  accordance  with  an  earlier  pro- 

29 

posal  by  Meyer -Peter  which  was  adopted  also  by  Taylor  and  Brooks  in  a 

30 

later  investigation.  Here,  S1  is  the  slope  that  would  be  required 
had  the  bed  been  flat  and  the  grain  size  fixed,  while  keeping  the  dis- 
charge the  same;  and  S2  is  the  additional  slope  as  a consequence  of 
the  bed  features  present.  Similarly 


Figure  21.  Coefficient  of  friction  of  plane  granular  surfaces  in  Fort  Collins  and  WES  tests 


Tb  T1  + T2 


and  as  before 


Tb  = PgRbSo’  T1  = PgRbSl’  T2  = pgRbS2 


The  resistance  coefficients  will  be  defined  as 


A A - 

Ab  y2  • A1 


2gVl 

u2 


> ^2 


2gRbS2 

u2 


(2.51) 


The  grain  resistance  coefficient  A is  preferably  determined  by  the 

1 27 

method  proposed  by  Lovera  and  Kennedy.  With  reference  to  A 2 , 
resort  is  made  first  to  dimensional  analysis.  Following  the  thesis  of 
Brooks,  S2  or  its  equivalent  A^  is  taken  as  a dependent  variable 


A2  = A2  (' 


U UR  o L 
r-r  ’ v ’ L ’ R ’ 


ps  \ 


(2.52) 


where 

k = grain  size  d^^ 

Ho  = height  of  dunes 
L = length  of  dunes 
Ps  = density  of  sand  grains 
p = density  of  water 

a = geometric  standard  deviation  of  sand  grains  in  the  bed 
The  authors  argue  that  the  resistance  of  ripples  and  dunes  is  not 
affected  by  viscosity.  Vanoni  and  Brooks  have  already  shown  that  \2 
is  not  affected  by  a change  in  a . The  ratio  Ps/P  is  constant  for 
ordinary  sand,  about  2.68.  Thus  Equation  2.52  may  simplify  to 


A2  ^2 


k o L 
* R * ~L  ’ R 


(2.53) 


In  accordance  with  the  Kennedy  theoretical  analysis  one  may  write 


■*(« 


(2.  5*+ ) 


V 


where  U is  the  critical  velocity  for  the  entrainment  of  sand  of  size 
C 2 

k . From  Shields'  diagram,  U#c/gk  is  constant  for  the  sands  commonly 
occurring  in  rivers.  Since  U/U#  is  a function  of  k/R  , this  required 


that 

and  thus  Equation  2.54  becomes 

H 


/gk 


(I) 


(2.55) 


o . / U k' 

L = *3  ’ R 


i>^k 


(2.56) 


Kennedy's  theoretical  analysis  suggests  also  that  the  length  of  dune  be 
given  as 


L = 2ir 


U 


g 


and  hence 


k - ow  uf.  k 

R ~ 2 gk  R 


(2.57) 


(2.58) 


Using  Equations  2.5**  and  2.57,  the  relation  for  the  resistance  coef- 
ficient from  Equation  2.53  is  simply 


A2 


(2.59) 


The  sufficiency  and  the  applicability  of  the  last  relation  was  examined 
by  Alam  and  Kennedy,  taking  a large  number  of  data  from  laboratory  chan- 
nels and  from  natural  rivers.  The  curves  of  Figure  22  summarize  the 
dependence  of  4x 2 = f2  on  the  relative  depth  R/k  for  varying  Froude 
numbers  U/ /gk  . These  curves  reveal  that  when  the  flow  Froude  number 
U//IR  is  smaller  than  0.5  the  resistance  coefficient  of  the  bed  fea- 
tures is  independent  of  R/k  and  is  essentially  a function  of  U//gk 

IJL 


(2.60) 


Alam  and  Kennedy  bring  out  the  point  that  the  parameter  ip'  of  the 
Einstein-Barbarossa  analysis  Equation  2.42  virtually  yields  the  relation 
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Figure  22.  Form  resistance  in  alluvial  channels  according 
to  Alam  and  Kennedy 


= 1-68  , k = d 

A1  u 


(2.C) 


As  X ^ does  not  vary  widely  among  different  alluvial  channels,  there 
is  not  much  variation  of  the  granulometric  distribution  of  the  natural 
sand  so  that  d^  is  proportional  to  d^^  . Accordingly  the  two  de- 
terminations of  X2  , one  by  the  Einstein-Barbarossa  method  and  the 
other  by  the  Alam-Kennedy  method,  would  be  close  to  each  other  for  the 
cases  where  Froude  numbers  are  less  than  0.5* 

55.  For  reference  it  may  be  desirable  to  have  an  expression  from 
the  Alam-Kennedy  study  for  X 2 when  Froude  number  is  less  than  0.5. 

The  averaged  values  of  X 2 over  the  horizontal  portions  of  the  curves 
in  Figure  22  are  plotted  against  U/ t^gk  in  Figure  23.  Representing 
the  alignment  of  the  points  by  a straight  line,  one  now  has 


79 


X2  = 0.5 


-3/2 


(2.62) 


If  we  introduce  the  relative  density  of  the  grains,  s = 2.68  , the 
expression  becomes 


A2  = 0.3U 


(2.63) 


Application  of  Alluvial  Channel  Resistance  Formulae 
to  Channels  Charged  with  Atlantic  Muck 


32 

56.  Kalinske  and  Hsia  performed  experiments  with  a bed  of  very 
finely  ground  silica  with  a specific  gravity  of  2.67.  The  gradation  of 
the  sand,  shown  in  Figure  2h,  is  similar  to  Atlantic  muck.  The  purpose 
of  the  experiments,  carried  out  in  a steel  channel  2.25  ft  wide, 

1.096  ft  deep,  and  80  ft  long,  was  to  study  the  transportation  of  such 
a fine  sediment  and  to  distinguish  between  bed  load  and  suspended  load. 
For  our  purpose  we  note  that  a range  of  flow  conditions  were  tested  and 
the  Manning's  n was  calculated  for  each  flow.  Data  associated  with 
ripples  in  the  bed  are  shown  in  Table  19-  The  entries  of  the  column 
for  Rq  , the  hydraulic  radius,  are  our  own  evaluations.  The  energy 
gradient  is  inferred  from  the  shear  velocity  given  in  the  last  column. 
We  shall  consider  this  data,  first  to  see  if  the  friction  of  the 
ripples  computed  by  the  Einstein-Barbarrosa  method  is  in  agreement  with 
the  curve  derived  by  these  investigators  for  the  American  rivers,  and 
secondly  to  inquire  if  the  values  obtained  are  in  agreement  with  the 
Alam-Kennedy  formula.  Equation  2.62. 

57.  Assuming  that  the  kinematic  viscosity  of  the  water  during 

-5  2 

the  Kalinske  and  Hsia  experiments  was  close  to  v = 1 x 10  ft  /sec 
and  the  smallest  value  of  the  friction  velocity  was  about  0.05  ft/sec, 
the  boundary  layer  6 would  then  equal  0.002U  ft.  Identify  k with 
dg,-  , that  is,  k = 0.012  mm  . The  ratio  of  the  boundary  layer  thick- 
ness to  the  grain  size,  6/k  , is  5^  • Hence,  we  are  dealing  with 
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3 

2.5 

0.37 

0.84 

0.278 

0.64 

0.0119 

0.046 

4 

2.5 

0.49 

1.14 

0.342 

I.29 

0.0106 

0.054 

5 

2.5 

0. 66 

1.31 

0.4l6 

1.67 

0.0105 

0.06l 

6 

5.0 

0.36 

I.25 

0.272 

1.95 

0.0117 

0.069 

T 

5.0 

0.52 

1.54 

0.354 

2.24 

0.0112 

0.080 

8 

5.0 

0.64 

1.8l 

0.410 

2.27 

0.0101 

0.084 

i 9 

10.0 

0.35 

1.68 

0.267 

3.36 

0.0119 

0.096 

1° 

10.0 

0.56 

2.14 

0.372 

6.81 

0.0114 

0.125 

hydrodynamically  smooth  surfaces  and  A to  be  inserted  in  Equation  2.38 
is  of  the  value 


A = 3.3 


U. 


*1 


The  resulting  equation  may  be  readily  transformed  to 


with 


= A + 8.65  log  (gRjS) 


A = 3.26  - 5.75  log  (gvS) 


(2.64) 


This  determines  the  hydraulic  radius  associated  with  the  surface 


and 


S , and  v are  known. 
R-^  is  known 

The  determinations  are 

*2 

(2.65) 

,'iV 

(2.66) 

i.67  d 

0,'  - 

R^S 

(2.67) 

58.  The  symbol  Rfc  appearing  in  Equation  2.65  denotes  the 
hydraulic  radius  of  the  channel  bottom  and  is  not  to  be  confused  with 
the  channel  hydraulic  radius  Rq  . In  channels  where  the  depth  H is 
negligible  in  comparison  with  the  channel  width  and  the  asperities  of 
the  sidewalls  are  smaller  than  those  of  the  bottom,  Rfe  would  equal 

Rq  ; otherwise,  Rb  would  be  greater  than  Rq  . The  matter  of  the  side- 
wall  corrections  was  discussed  before.  As 


from  Equation  2.24  in  a rectangular  channel 
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•v,*~  V 


2 -1/3 

and  in  view  of  Equation  1.20  and  as  A ~ n R 


/n\3/2 

v-k) 


where  n^  is  the  Manning's  n of  the  roughness  of  the  sidewalls,  n , 
the  Manning's  n of  the  channel 


l 


(2.68) 


Kalinske  and  Hsia  note  that  n of  the  still  channel  is  0.0092.  The 

w 

Manning's  n's  of  the  channel  with  the  ripples  present  are  shown  in 
Table  19-  Determinations  of  the  procedures  in  Equations  2.64  to  2.66 
are  shown  in  Table  20.  Values  of  the  parameters  appearing  in  the  last 
two  columns  of  the  table  are  plotted  in  Figure  l8.  Ignoring  the  scatter. 


Table  20 


Determination  of  the  Einstein  and  Barbarossa  Variables 


in  the  Kalinske  and 

Hsia  Data 

(Assumed: 

v_ 

= 1 x io“ 5 

ft^/sec ) 

Run 

No. 

4 

S X 10 

U , ft/ sec 

«b  • ft 

Rx  , ft 

r2  , ft 

u/u*2 

3 

2.36 

0.84 

0.306 

0.226 

0.080 

34.0 

0.77 

4 

2.64 

1.14 

0.369 

0.323 

0.046 

58.0 

0.48 

5 

2.97 

1.31 

0.461 

0.385 

0.076 

50.2 

0.38 

6 

5.56 

1.25 

0.300 

0.199 

0.101 

29.3 

0.37 

7 

5.61 

1.54 

0.386 

0.267 

0.119 

33.2 

0.27 

8 

5.30 

1.81 

0.443 

0.361 

0.082 

48.4 

0.21 

9 

11.98 

1.68 

0.272 

0.201 

0.071 

32.2 

0.17 

10 

13.00 

2.14 

0.421 

0.218 

0.203 

26. 6 

0.15 

Note:  S = surface  energy  gradient  computed  using  the  entries  of  the 
last  column  in  Table  19,  T0  = PS^S  . 


84 


« <' 

* n 


t 


the  average  values  fall  on  the  curve  from  the  Einstein-Barbarossa 
equation 


yp-  = 30UT1/2 
u#2 


(2.1*3  bis) 


The  scatter  of  the  plotted  points  is  very  severe.  Probably  this  may  be 
due  to  the  difficulties  in  measuring  the  energy  gradient  of  the  flows. 
Computations  made  on  the  basis  of  the  experimental  channel  slope  show  a 
lesser  degree  of  scattering.  The  method  is  also  very  sensitive  to 
slight  inaccuracies  in  the  determination  of  the  current  mean  velocities. 

59.  Since  all  the  Froude  numbers  of  the  flows  in  the  Kalinske 
and  Hsia  tests  are  less  than  0.5,  the  coefficient  X ^ of  the  features, 
consisting  essentially  of  ripples,  may  be  also  obtained  by  the  formula 


= 0.5^—  ) 

Vik  / 


(2.62  bis) 


where  k , the  median  grain  size,  has  the  value  of  0.010  mm.  It  will  be 
recalled  that  the  expression  was  based  on  the  Alam-Kennedy  analysis 
curves  giving  the  dependence  of  X 2 on  H/k  and  U/gk  . The  determi- 


nations of 


X2  on  this  basis  is  shown  in  Table  21.  In  the  same  table 


are  given  also  the  corresponding  values  in  accordance  with  the 
Einstein-Barbarossa  method.  The  formula  used  was 


' 2(u*2) 


and  the  appropriate  values  were  taken  from  the  seventh  column  in 
Table  20.  Although  the  agreement  of  the  values  for  the  individual  runs 
determined  by  the  two  methods  is  not  always  satisfactory,  the  agreement 
of  the  means  is  good.  We  now  conclude  that  either  method  would  be 
satisfactory  for  the  determination  of  X of  the  bed  features  consist- 
ing of  ripples  or  dunes  of  flows  with  a Froude  number  less  than  0.5. 

60.  The  determination  of  the  resistance  of  a deep  and  wide  chan- 
nel charged  with  Atlantic  muck  may  be  made  using  the  principles  dis- 
cussed in  the  above  sections.  At  this  point,  however,  it  would  be 
useful  to  remark  that  the  Einstein-Barbarossa  method  would  be  equally 
applicable  if  the  surface  energy  gradient  S is  broken  into  two  parts. 


• * * mi  VJ*V*'*’  V ■ 
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, and  Sg  , the  first  associated  with  the  surface  friction  and  the 
second  with  the  form  resistance  of  the  bed  features  in  the  form  of  rip- 
ples and  dunes.  We  note  that  R = 

RXS  = RS^ 


and 


R2S  = R1S2 


considering  a hydrodynamically  smooth  surface  from  Equations  2.38  and 

2.45 

U, 


U = /T7 

u*i  "V  - 


= 3.25  + 5.75  log  -jp  + 5.75  log  ^ (2.38  bis) 


where 


U#1  = 


Thus  U/U*.  is  a function  of  UR/v  and  the  relation  of  the  function 
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was  shown  in  Table  17.  After  specifying  the  values  of  , R , 

and  v the  value  of  is  determined  by  means  of  the  above  relation 

The  gradient  S0  is  obtained  from 


now 


Reductions  from  laboratory  and  alluvial  river  data  made  on  this  basis 
would  not  change  the  results  already  shown  in  Figure  l8  since  RS  = ! 


For  the  present,  however,  it  would  be  easier  to  revert  to  Alam-Kennedy 

procedure  to  determine  the  form  resistance  of  a deep  channel.  In  the 

channel  previously  considered  in  paragraph  31,  R = 50.6  ft  and 
-5  2 

v = 1 x 10  ft  /sec.  Now  n can  be  computed  from 


after  determining  A from  Equation  2.38  or  using  Table  17.  The  re' 
suiting  values  agree  with  n given  in  paragraph  31.  As  regards  n0 


and  A_  may  be  computed  from 


>2 = °-5fe) 


(2.62  bis) 


for  U/f'gR  less  than  0.5.  For  Atlantic  muck  k = 0.005  mm.  The 
combined  n for  this  channel  is  as  follows: 


ft/ sec 

nl 

n2 

n 

0.5 

0.0159 

0.0262 

0.0421 

1.5 

0.0152 

0.0146 

0.0298 

2.0 

0.0146 

0.0088 

0.0234 

4.0 

0.0139 

0.0052 

0.0191 

6.0 

0.0136 

0.0039 

0.0175 

8.0 

0.0131 

0.0031 

0.0162 

10.0 

0.0130 

0.0026 

0.0152 

The  Manning's  n of  the  channel  considered  decreases  when  the  current 
velocity  is  increased.  In  the  above,  A^  was  calculated  on  the  assump- 
tion that  the  surface  friction  is  that  of  a hydrodynamically  smooth 
surface.  If  we  had  used  the  method  of  Lovera-Kennedy  for  establishing 
A^  , these  would  have  recalled  n values  larger  than  those  shown  in 
the  above.  Judging  from  the  curves  shown  in  Figure  21,  for  moving 
particles  A^  is  almost  twice  as  large  as  if  the  particles  had  been 
fixed.  This  would  increase  the  above  shown  values  of  n by  l.1*  times. 


PART  III:  ANALYSIS  OF  FORM  RESISTANCE  IN  LABORATORY  CHANNELS 


6l.  A second  method  for  obtaining  the  resistance  of  alluvial 
channels  would  be  to  relate  the  form  resistance  of  the  bed  features 
present  to  the  dimensions  of  the  latter.  Probably  this  will  be  pos- 


sible more  readily  to  features  consisting  of  ripples  and  dunes.  In 
principle,  then,  a formula  needs  to  be  developed  relating  the  resis- 
tance coefficient  A 2 to  the  height  and  length  of  the  dimes  or  the 
ripples.  In  addition,  formulae  need  to  be  developed  relating  the  size 
and  the  length  of  the  features  to  grain  size  and  channel  stress  or 
channel  current  velocity. 

62.  Let  F^  be  the  resistive  force  of  the  ripples  or  dimes 

per  unit  width  of  channel  and  H the  height  of  the  dimes  or  ripples. 

8 ® 

Using  the  formula  of  Sadron 


F 

r 


CDHo  t U 


2 

o 


where 

CD  = the  coefficient  of  drag  force 

p = the  density  of  the  flowing  liquid 

U = the  current  at  the  crests  of  the  bed  features 
o 

Introducing  x2  the  stress  arising  from  the  bed  features  and  L the 
length  of  the  ripples  on  the  dunes 


F 


r 


o H 

c u2  — 

Do  2 L 


Dividing  by 


, the  channel  current  velocity  squared 


then  expressing  the  water  depth  by  H 


Substituting  this  in  Equation  3»1  and  arranging  terms 


Thus  C_  may  be  established  from  observed  values  of  H H 

u o’’ 

L provided  m , the  exponent  in  the  velocity  law  (Equation  3.2) 
known. 


Observed  velocities 
in  Fort  Collins  tests 


63.  In  the  Fort  Collins  tests 
approximated  by  the  form 


where  C is  the  velocity  at  a point  1 ft  from  the  bottom  of  flume 
This  expression  can  be  changed  to  power  law 


Comparing  these  expressions  with  each  other 


This  is  satisfied  for  y 
y = 0.1  ft.  Accordingly 


Assume  that  it  is  satisfied  also  for 


I 


m = -log  (l  - §) 

The  constants  C and  B are  given  in  the  appropriate  tables  of  ref- 
erence. Through  the  latter  relation,  n was  evaluated  for  ripples 
and  also  for  dunes.  It  was  noted  that  m did  not  show  much  variation 
from  test  to  test.  Averaged  values  are  shown  in  Table  22,  where  the 

Table  22 


Exponent  m of  Velocity  Distribution  in  Fort  Co..lins  Data 

(u  = Cym) 


Grain  Size 
d50 

Ripple 

m 

Dune 

m 

Average 

m 

0.19  mm 

0.196 

0.151 

0.173 

0.27  and  28 

0.266 

0.116 

0.191 

0.1*5 

0.161* 

0.133 

0.11+8 

0.1*7 

0.183 

0.158 

0.170 

0.93 

— 

0.137 

0.137 

Average 

0.161* 

entries  are 

distinguished  both  as  grain  size  and  feature  forms.  The 

exponent  m 

is  practically  independent  of 

grain  size.  Its 

value  is 

somewhat  greater  for  the  ripples.  Ignoring  this  difference  we  adopt 
the  average 


m = 0.166 


1 


: 


As  an  example,  the  distribution  of  velocities  over  a ripple  and  also 
over  a dune  is  shown  in  Figure  25.  In  a cross  section,  the  velocities 
were  measured  at  the  midpoint  and  at  points  close  to  the  two  walls. 
Using  this  value  of  m.  Equation  3.3  becomes 


T 


T 


T 


Figure  25.  Velocity  distribution  over  a ripple  and  a dune 
(8-ft  channel,  d^Q  = 0.19  mm) 

Determination  of 
for  Fort  Collins  data 

64.  The  energy  gradient  S is  divided  into  two  parts,  S 
balancing  the  resistance  due  to  granular  bottom  surface  and  S2 
ancing  the  form  resistance  of  the  ripples  or  of  the  dunes.  Thus 


Tb  = T1  + T2 


with  the  usual  definitions 


Tb  = pesV  Ti  = p«siV  t2  = Pgs^ 


U*  = gsi^;  = gs^;  U^2  = gS^ 


and  now  since  S = S + 

2 1 ^ 

Dividing  by  U and  rearranging  terms 
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(u*2)‘ 


(U,)‘ 


<V‘ 


of  the  terms  on  the  right-hand  side  the  first  is  taken  from  the  observe- 

n Q 

tion  of  the  Fort  Collins  data  after  sidewall  corrections.  As  regards 
the  second,  this  was  determined  in  two  ways  using  first,  the  ordinary 
asperities  formulas,  putting  H equal  to 


C 5.751CS  12.8 
*1  s 


where  kg  stands  for  the  grain  median  diameter 


u , U*1H 

rr-  = -3.^  + 5. 75  log  — ■ — 

v 


(2.38  bis) 
and  secondly, 

(2.50  bis) 


It  will  be  recalled  referring  to  Figure  21  that  the  latter  expression 

represents  on  the  average  the  granular  skin  friction  found  from  the  Fort 
l8  2 

Collins  tests.  After  evaluating  (U#2)  /U  for  each  run  and  noting  the 
quantities  Hq/H  and  Hq/L  , these  are  entered  into  Equation  3.1*  to 
obtain  C^  . The  computed  C^'s  by  these  two  methods  are  plotted 
against  L/Hq  in  Figures  26  and  27  in  order  to  find  if  there  is  de- 
pendence on  the  relative  lengths  of  the  bed  features.  The  data  are 
differentiated  ps  to  the  grain  size  and  bed  feature  form.  No  depen- 
dence on  the  grain  size  is  indicated.  The  ordinates  represent  the 
averages  of  C^  corresponding  to  L/Hq  constant  or  very  close  to  each 
other.  There  is  much  scatter  which  is  expected  owing  to  the  difficul- 
ties in  measuring  H and  the  fact  that  represents  a difference 

O ^d 

of  terms,  each  open  to  error.  It  is  also  possible  that  L and  Hq 

alone  by  themselves  are  not  sufficient  to  formulate  the  stresses  of  the 

33 

bed  features.  This  problem  has  been  discussed  by  Nordin.  Owing  to 
the  scatter  of  the  points,  it  is  difficult  to  give  an  exact  graphics  to 
define  the  dependence  of  C^  on  L/Hq  . Tentatively,  the  horizontal 
lin^s  are  drawn  in  the  two  figures  implying  that  C^  is  independent  of 
L/Hq  . From  Figure  26,  which  is  the  case  in  which  the  granular  surface 
resistance  is  computed  by  Equation  2.38,  the  roughness  formula. 
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Figure  26.  Dependence  of  the  coefficient  of  form  resistance  on  feature 
sizes  and  water  depth  in  Fort  Collins  8- ft  flume  tests.  Granular 
surface  friction  computed  by  ordinary  formula  (Equation  2.38) 


and  introducing  into  Equation  3.4  we  find 


From  Figure  26,  which  is  the  case  in  which  the  granular  surface 
resistance  is  computed  by  Equation  2.50,  modified  Prandtl  formula 


Figure  27.  Dependence  of  the  coefficient  of  form  resistance  on  feature 
sizes  and  water  depth  in  Fort  Collins  8- ft  flume  tests.  Granular 
surface  friction  computed  by  modified  Prandtl  formula  (Equation  2.50) 


and  introducing  into  Equation  3. ^ we  find 


and 


(3.7) 


Here  in  Equations  3.6  and  3.7  we  have  two  different  expressions  for  the 

resistance  coefficients  of  bed  features  in  terms  of  H /H  and  H /L 

o o 

which  illustrate  the  importance  of  the  law  of  the  granular  surface 

friction  for  the  determination  of  the  coefficient. 

65.  The  above  conclusion  that  is  independent  of  L/ Hq  in 

the  Fort  Collins  data  is  based  on  the  premise  of  (U/U  )^  being  pro- 
l/3  ® 

portional  to  (Hq/H)  ' . Conversely,  it  would  be  important  to  find  if 


the  assumption  of  Cn  being  independent  of  L/H  would  lead  necessarily 

u 2 ® -\  /o 

to  the  relation  that  (U/U  ) is  proportional  to  (H  /H  ) ' j . Since 

o o o 

CD  is  constant,  we  write  Equation  3.3  in  the  form 


(3.8) 

(3-9) 


Computing  of  the  Fort  Collins  data  on  the  basis  of  Equation  3.9,  the 

results  are  plotted  against  Hq/H  in  Figures  28  and  29.  The  first 


Figure  28.  Dependence  of  the  coefficient  of  form  resistance  on  water 
depth  in  Fort  Collins  8-ft  flume  tests.  Granular  surface  friction  com- 
puted by  roughness  formula  (Equation  2.38) 

represents  the  case  in  which  the  granular  surface  friction  is  computed 
by  Equation  2.38  and  the  second  where  the  granular  surface  friction  is 
computed  by  Equation  2.50.  The  data  are  differentiated  as  to  the  grain 
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Figure  29.  The  dependence  of  the  coefficient  of  form  resistance  on  water 
depth  in  Fort  Collins  8-ft  flume  tests.  Granular  surface  friction  com- 
puted by  modified  Prandtl  formula  (Equation  2.50) 


size  and  bed  feature  form.  No  dependence  on  grain  size  is  indicated. 
Points  represent  the  average  value  in  initial  data  of  Hq/H 

being  constant  or  close  to  each  other.  Although  there  is  scatter  of 
points,  the  scatter  is  not  as  severe  as  that  in  Figures  26  and  27.  The 
lines  of  the  figures  are  drawn  such  that  the  number  of  the  points  lying 
above  the  line  is  the  same  as  those  lying  below.  The  lines  yield 


and 


and  from  Equation  3.8  we  find  that  the  corresponding  values  of  C_^  are 


0.44  and  0.34.  These  are  in  agreement  with  the  ordinates  of  the  lines 
shown  in  Figures  2 6 and  27. 

66.  In  Figures  26  and  27,  the  average  value  of  for  the 

dunes  is  larger  than  the  average  value  of  for  the  ripples.  This 

suggests  the  possibility  that  although  CD  is  independent  of  L/Hq  , 

it  has  a lesser  value  for  the  ripples  than  for  the  dunes.  To  examine 

the  point  we  give  new  values  for  C from  the  ripple  data  of  Vanoni  and 

34  D 

Hwang  in  Figures  30  and  31  and  also  some  results  from  WES  (Franco 
tests).  The  latter  results  have  not  been  published  and  are  shown  in 


7 0.23  FRANCO 

Q 0.23  VANONI 
O 0.14  VANONI 
• 0.21  VANONI 


Figure  30.  Dependence  of  the  coefficient  of  form  resistance  on  ripple 
size  and  water  depth.  Granular  surface  friction  computed  by  ordinary 

formula  (Equation  2.38) 

Table  23  with  the  kind  permission  of  Mr.  Franco.  The  geometry  of  the 
ripples  had  been  examined  for  a possible  dependence  of  the  size  on 
water  temperature.  The  surface  displacement  of  the  bottom  sand  had 
been  surveyed  from  one  end  of  the  flume  to  the  other.  The  height  and 
the  length  of  the  ripples  were  deduced  from  the  extremes  in  the 
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Figure  31.  Dependence  of  the  coefficient  of  form  resistance  on  ripple 
size  and  water  depth.  Granular  surface  friction  computed  hy  modified 

Prandtl  formula  (Equation  2.50) 

variation  and  the  mean.  Horizontal  lines  are  drawn  in  both  Figures  30 

and  31.  The  suggestion  is  that  CL  is  independent  of  L/H  . The 

u o 

lines  are  drawn  such  that  the  number  of  points  above  the  line  is  equal 
to  the  number  of  points  below  the  line.  These  give 


CD  = 0.23 

if  the  surface  granular  friction  is  computed  by  Equation  2.38  roughness 
formula,  and 

CD  = 0.18 

if  computed  by  Equation  2.50,  modified  Prandtl  formula.  Accordingly, 
in  place  of  Equation  3.6 


a ^ i?  art  i 
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Table  23 


Ripple  Dimensions  of  WES  (Franco  Tests ) 


U 

U*b 

*b 

V 

p 

H 

0 

L 

*s 

No. 

ft/ sec 

ft/sec 

ft 

ft  /sec 

ft 

ft 

H = 0.5  ft 

21 

1.197 

0.1363 

0.1+83 

1.67  x io-5 

0.0756 

0.6435 

5.6 

h 

1.183 

0.1319 

0.1+87 

1.67 

0.0754 

0.5913 

lb/hr/ft 

6 

1.192 

0.1301+ 

0.1+85 

1.21 

0.061+6 

0.5078 

5 

1.198 

0.1201+ 

0.1+80 

0.93 

0.0632 

0.4280 

15 

1.294 

0.11+71 

0. 1+1+9 

1.67 

0.0832 

0.6915 

13.3 

10 

1.312 

0.1370 

0. 1+1+2 

1.21 

0.071+1+ 

0.6724 

14 

1.352 

0.1268 

0.1+98 

0.93 

0.0710 

0.7581 

I 

1 

11 

1.368 

0.1523 

0.1+25 

1.67 

0.0888 

0.6373 

23.3 

12 

1.1+09 

0. 11+22 

0.1+11 

1.21 

0.0776 

0.7547 

S 

10 

1.1+22 

0.1278 

0.1+03 

0.93 

0.0772 

0.6863 

!!  ' 

17 

1.1+76 

0.1278 

0.388 

0.93 

0. 0726 

0.8064 

* 

H = 0.75  ft 

1 

30 

1.513 

0.1526 

0.51+0 

1.67  x 10-5 

0.07I+8 

0.6309 

23.3 

• - 

| 4 

31 

1.61+0 

0.1332 

0.1+89 

0.93 

0.0706 

0.6644 

■ 

! 

32 

1.31+5 

0.1331 

0.598 

1.67 

0. 0684 

0. 6031 

5.6 

n 

3b 

1.31+5 

0.1198 

0. 59U 

0.93 

0.0644 

0.5530 

r y 

[ h 

33 

1.354 

0.121+8 

0.591 

1.21 

0.0652 

0.5518 

\\ 

H (l 

a/3 

I 

a2  = 0. 

" f2 

(3.10) 

! 

I 

and 

in  place  of  Equation  3. 7 

! > 

H /} 

a/3 

r 

! 

a2  = 0. 

2hJt\- 

i)  S 

= f2 

(3.11) 

These  apply  to  ripples.  Interestingly  enough,  the  last  expansion  is 
very  close  to  the  Kronoz  formula 

0.25 

(3.12) 


H /H  \ ' 

x2  = 0*21  “L  vi) 


cited  by  Jensen  and  Lebreton. 
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Dependence  of  Ripple  and  Dune  Size  on 
Current  Mean  Velocity 


67.  It  is  natural,  intrigued  by  the  Alam  and  Kennedy  simple 

results  on  the  resistance  of  alluvial  streams,  to  inquire  if  similar 
results  could  be  obtained  applicable  also  to  laboratory  flumes.  Accord- 
ing to  Alam  and  Kennedy,  when  the  Froude  number  is  less  than  0. 5 the 
form  resistance  coefficient  in  alluvial  rivers  depends  solely  on 

U//gk  where  k is  the  median  sand  diameter  d^  . For  the  desired 

18 

examination  we  shall  make  use  of  that  part  of  the  Fort  Collins  data 
that  relates  to  ripples  and  dunes.  These  were  observed  for  flows  with  a 
Froude  number  less  than  0.5-  In  the  preceding  analysis  X 2 is  ex- 
pressed in  terms  of  Hq/L  and  Hq/H  . Hence,  it  is  required  that  the 

dependence  of  these  on  U2/(s  - l)gd  be  established. 

2 

68.  The  dependence  of  HQ/d  on  U /(s  - l)gd  in  the  Fort  Col- 
lins data  is  indicated  in  Figure  32.  The  points  are  differentiated  as 
to  grain  size  and  bed  feature  form.  The  ordinates  in  some  cases 

represent  the  average  for  the  points  having  a constant  or  close  abscissa. 

2 

The  dependence  of  HQ/d  on  U / ( s - l)gd  for  the  dunes  differs  consid- 
erably from  that  of  the  ripples.  For  the  dunes 

-j  = 1*0  [u2/(s  - l)gd]0,1*5  (3.13) 

and  for  the  ripples 


-j-  = 10.3[u2/(s  - l)gd]°*45  ( 3.  i*0 

if  some  of  the  points  of  the  tests  with  0. U-mm  sand  are  disregarded. 

For  the  dunes  these  lead  to  the  expression 

-2-=  l+0[u2/(s  - l)gd]0,1*5  | (3.15) 


and  for  the  ripples 


Figure  32.  Dependence  of  bed  feature  heights  on  current 
velocity  in  Fort  Collins  tests  (8- ft  flume) 

= 10.3  [u2/(s  - l)gd]0‘U5  | (3.16) 


The  last  relation  furnishes  the  explanation  why  in  the  models  of  water- 
courses the  ripple  sizes  are  exaggerated.  The  relation  can  be  written  as 


10.3  [u2/(s  - l)gH]1/2(|) 


(3.17) 


The  model  currents  being  determined  by  Froudian  relation,  the  ratio 
of  the  ripple  heights  in  the  model  and  in  the  prototype  are 


H 


op 


(3.18) 


69.  Next  the  dependence  of  L/Hq 
Collins  data  is  indicated  in  Figure  33. 


2 

on  U / ( s - 1 ) gd  in  the  Fort 
The  points  are  differentiated 
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Figure  33.  Dependence  of  the  length-to-depth  ratio  of  bed  features 
on  current  velocity  in  Fort  Collins  tests  (8- ft  flume) 

as  to  grain  size  and  the  bed  feature  form.  In  some  cases,  the  ordinate 
values  represent  the  averages  for  the  points  having  a constant  or  close 
abscissa.  Apparently,  the  dependence  of  L/Hq  on  U /(s  - l)gd  is  the 
same  for  the  ripples  and  for  the  dunes.  The  relation  seems  to  be  prac- 
tically independent  of  the  grain  size.  Accordingly 


(3.19) 


if  some  of  the  points  of  the  tests  with  0.1+5-mm  sand  are  disregarded. 

70.  The  substitution  of  H /H  from  Equation  3.15  or  3.16  and  of 
L/H  from  Equation  3.19  in  Equation  3.6  or  3.7  should  yield  the  ex- 
pressions of  the  form  resistance  coefficient  Xg  for  dunes.  Thus  for 
the  dunes,  putting  S = 2.68  , 


X_  = 1.80(U//gd)  "1,3°  (d/H)1/3  , X0  = kt0 


(3.20) 


using  Equation  3.6  (granular  friction  computed  by  roughness  formula),  on 


A2  = 1.28(U/i^gd)  1‘3°(d/H)l/3 


(3.21) 


using  Equation  3.7  (granular  friction  computed  by  modified  Prandtl 
formula).  As  regards  the  ripples 


X2  = 0.79(U//id)_1'3O(d/H)1//3 


(3.22) 


using  Equation  3.10  (granular  friction  computed  by  ordinary  roughness 
formula ) , or 


X2  = 0.6l(U/v^d)  1-30(d/H)l/3 


(3.23) 


using  Equation  3.11  (granular  friction  computed  by  modified  Prandtl 
formula).  Now  there  is  this  puzzling  difference  in  the  behavior  of 
laboratory  flumes  and  the  natural  watercourses.  While  in  the  natural 
rivers  the  resistance  coefficient  is  practically  independent  of  the 

ratio  d/H  and  depends  solely  on  U/i/gd  for  flows  with  a Froude  number 
less  than  0.5,  in  flumes  there  is  also  dependence  on  the  ratio  of  d/H 
even  though  the  Froude  number  of  currents  is  below  0.5.  The  graphs  of 
Equation  3.20  for  the  H/d  = 1 * 103  and  H/d  = 0. 5 x 103  are  shown 
in  Figure  3^  and  the  points  shown  between  the  two  curves  are  for  rivers 
according  to  the  Alam  and  Kennedy  analysis. 

71.  The  expression  of  L/Ho  by  Equation  3.19  can  be  obtained  in 
another  way.  Values  of  the  ratio  L/d  computed  from  the  Fort  Collins 

1 Q 

data  are  shown  in  Figure  35.  The  data  are  identified  as  to  grain  size 
and  bed  feature  form.  The  ordinates  of  some  represent  the  average  for 
the  points  having  constant  or  very  close  abscissa  values.  The  alignments 
of  the  points  for  the  dunes  and  the  ripples  differ  from  each  other.  For 


the  dimes 


= U5.0  [u2/(s  - l)gd]5A 


(3.2k) 
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Figure  34.  Comparison  of  form  resistance  coefficient  based  on 
Fort  Collins  tests,  8- ft  flume,  with  Alam-Kennedy  river  values 


and  for  the  ripples 


| = 12.1+  [u2/(s  - l)gd] 


if  some  of  the  points  of  the  tests  with  the  O.i+5-mm  sand  are  disre- 
garded. These  equations,  in  conjunction  with  Equations  3.13  and  3-1^ 
lead  nearly  to  Equation  3.19. 


Dependence  of  bed  feature  lengths  on  current  velocity- 
in  Fort  Collins  tests  (8- ft  flume) 


72.  The  literature  dealing  with  sediment  problems  was  examined  to 

see  if  there  have  been  studies  that  have  dealt  with  the  question  of 

dependence  of  the  ratio  of  feature  lengths  to  the  size  of  grains  on 

Froude  number.  Most  combinations  on  the  geometry  of  bed  features  have 

35 

ignored  this  question,  except  one  given  by  Jensen  and  Lebreton.  Ac- 
cording to  these  authors,  based  on  tests  with  sediments  of  varying 


densities  and  with  coarse  materials  satisfying  the  condition  d#  > l6  , 
the  relation  is 

j = 193 (s  - l)0,1*8(U//gH) 

This  can  be  written  after  replacing  0.48  with  0.50  as 

£ = 193  [u2/(s  - l)gH]"l/2 

or 

193  [u2/(s  - l)gd]~1/2  (d/H)1/2  (3.26) 

a result  which  is  so  markedly  different  from  the  relation  expressed  by 
Equation  3.24.  In  the  Fort  Collins  experiments,  the  L/d  of  dunes  in- 
creases with  increasing  Froude  numbers.  In  great  contrast  with  this, 
in  the  Chatou  tests  L/d  decreases  with  increasing  Froude  number.  As 
in  the  Fort  Collins  tests,  d^  is  smaller  than  16  and  in  the  Chatou 
tests  it  is  larger  than  l6;  then  the  question  arises  if  these  conditions 
do  have  a bearing  on  the  above  difference  noted. 

A Discussion  on  the  Formulas  of  Form  Resistance 
and  Feature  Geometry 

73.  We  certainly  wish  to  consider  a reputed  formula  of  form 

34 

resistance  for  ripples  given  by  Vanoni  and  Hwang 

-==  - 3.5  log  ^ - 2.3  (3.27) 

/4A2  o 

where 

= hydraulic  radius  of  the  bed 
e = exposure  ratio  AgA 

A^  = horizontal  projection  of  the  lee  forces  of  the  ripples 

A = total  area  under  which  A is  measured 

s 

Hq  = mean  height  of  the  ripples 

As  the  heights  of  the  ripples  are  controlled  by  the  angle  of  repose  of 
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the  grain  particles,  A is  proportional  to 

s 

the  crests  of  ripples.  One  may  then  take  e 
and  write 


e = 


L , the  spacing  between 
to  be  proportional  to  Hq/L 


In  the  data  considered  by  Vanoni  and  Hwang,  a is  close  to  unity  and 
subject  to  variation  from  one  flume  to  another.  In  the  data  of  a bO- ft 
channel,  the  mean  value  of  a is  1.11;  of  a 60-ft  channel,  O.96;  and 
of  a 130-ft  channel,  I.36.  The  average  is  a = l.lU.  Again  in  the  same 
ratio  H/R^  also  varies  slightly  from  flume  to  flume  with  an  average 
of  H/R^  =1.09  . Inserting  these  quantities,  the  resistance  formula 
reduces  to 


— — = 3.5  log  ^ - 2.0 
AIT  H_ 


This  may  be  expressed  as  a power  law 


\2  = 0.1*15 


(3.28) 


(3.29) 


valid  for  the  range  of  Hq/L  x Hq/H  from  20  to  150,  the  range  covered 
by  the  data  considered  by  Vanoni  and  Hwang.  In  arriving  at  the  formula, 
these  authors  have  computed  the  granular  skin  friction  coefficient  A 
using  the  ordinary  roughness  laws.  Therefore,  the  above  expression 
should  be  compared  with  A 2 given  by  Equation  3.10;  that  is 

H /H  \0-33 

A2  = 0.31  (3.10  bis) 

a formula  of  the  similar  form  but  with  different  exponents  of  the  cor- 
responding quantities.  The  larger  difference  is  in  the  exponent  of 
Hq/H  . It  will  be  recalled  that  the  latter  was  developed  assuming  that 

the  velocity  over  the  forms  varies  as  y with  m = 0.166  in  the  data 

in 

of  Fort  Collins  tests.  Equation  3.29  suggests  that  for  the  ripples  con- 
sidered by  Vanoni-Hwang,  m can  be  as  high  as  0.38.  To  examine  this 
matter,  the  dependence  of  the  quantity  A2/(Hq/L)  on  Hq/H  both  for 
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the  Franco  data  and  also  of  the  data  considered  by  Vanoni-Hwang , Fig- 
ure 36  was  prepared.  Individual  points  represent  averages  from  the 
original  data  with  close  abscissa  values  taking  two  or  three  at  a time. 
Despite  the  scatter,  it  is  possible  to  represent  the  alignment  of  values 
from  the  Vanoni-Hwang  data  by  a line  drawn  in  the  figure  leading  to 


VANONI  d 0.23  mm  40-FT  FLUME 
VANONI  d 0.0137  mm  60-FT  FLUME 
VANONI  d 0.206  mm  13O.FT  FLUME 
FRANCO  d - 0.23  mm  H=  0.5  FT 
FRANCO  d = 0.23  mm  0.7  FT 


0.4  0.6  0.8  1 


Figure  36.  Dependence  of  coefficient  of  form  resistance  of 
ripples  on  feature  dimensions 


H /H  \ ‘ 

= 0^-r(i) 


(3.30) 


This  is,  in  a way,  closer  to  Equation  3.10  and  suggests  that  for  the 
ripples  considered,  m amounts  to  0.23.  The  data  from  Franco  tests  in- 
dicate a parallel  alignment.  The  implication  of  the  present  results  is 
that  there  is  a difference  in  the  form  resistances  of  ripples  and  dimes, 
a matter  that  we  shall  be  considering  on  another  occasion. 

74.  Jensen  and  Lebreton  state  that  the  Vanoni  formula  applies 
equally  well  to  their  experimental  results.  The  formula  they  give  is 
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(3.31) 


. 3.83.0,  (|-.y.2.76 

which  is  based  on  experiments  carried  out  in  an  inclined  channel  50  cm 

deep,  80  cm  wide,  and  21  m long.  The  grains  used  were  of  different 

densities  and  sizes,  the  latter  varying  from  1.2  to  2.4  mm.  The  nondi~ 

mensional  diameter  d#  exceeds  the  critical  value  16.5-  Here,  no 

2 

ripples  could  be  generated  and  as  the  stress  parameter  F#  was  greater 
than  0.05,  dunes  were  formed.  In  the  tests,  on  the  average  R^  = 0.88H 
and  accordingly  the  last  relation  reduces  to 

— — « 3.83  log(§“  • §-)-  2.55  f2  = U\2  (3.32) 

Al2  V o o / 

Expressing  this  as  a power  formula 

/H  \°-kk/I,?-kh 

*2  - 0.076  (-^j  yj  f2  = U2  (3.33) 


2 —1 

valid  in  the  range  (Hq/LH)  from  100  to  2000,  the  range  covered  by 
the  Jensen  and  Lebreton  data.  These  authors  give  an  additional  ex- 
pression in  the  form 


(3.34) 


The  constants  , n , and  m are  determined  by  the  method  of  least 
squares.  Results  of  the  Jensen-Lebreton  data  may  be  changed  to 


Is  \°-66 

X2  - 0.107\-f;  \±)  f2  = U2  (3.35) 

Equations  3.32  and  3.35  as  applying  to  the  same  experimental  data  should 
have  been  identical,  having  the  same  factor  A and  the  same  exponents 
m and  n . Owing  to  the  scatter  in  the  data,  this  condition  is  not 
realized. 


-t*  * v«* 
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75.  A criticism  generally  expressed  regarding  the  usual  formula 
of  resistance  of  alluvial  streams  and  of  the  geometry  of  features  is  that 
factors  expressing  the  influence  of  viscosity  are  lacking.  The  criticism 
is  on  the  grounds  that  a change  in  the  temperature  of  waters  affects  the 
resistance  of  alluvial  streams.  The  problem  of  the  influence  of  the 
effect  of  viscosity  is  well  documented  by  Vanoni"5  by  examining  the  im- 
portance of  the  i/gd  * d/v  for  the  prediction  of  the  bed  features. 

37 

Simons  and  others  find  that  the  Barekyan  formula  is  quite  applicable 
to  the  determination  of  dune  heights  in  the  Fort  Collins  data,  yet  in 
this  formula  an  appropriate  viscosity  parameter  is  missing.  In  the 
presentation  we  now  give  in  Figure  37  this  defect  is  remedied  by 
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Figure  37.  Ripple  height  relative  to  grain  median  diameter  in 
Fort  Collins  and  WES  tests 


introducing  the  grain  parameter  d#  which  involves  the  kinematic  vis- 
cosity v . In  the  figure,  H /(dd„)  is  plotted  against  where 

o w o 

i|<o  is  the  Shields'  parameter  for  the  initiation  of  sand  entrainment. 
Included  also  in  the  figure  are  points  from  the  WES-Franco  data  which 
are  in  good  agreement  with  the  Fort  Collins  data.  On  the  present  basis, 
the  expression  for  the  ripple  height  is 

HQ  = U5(*  - *0)dd*  (3.36) 


where 


76.  Owing  to  the  similarity  in  the  scatter  of  data  in  the  two 
representations,  one  in  accordance  with  the  Barekyan  plan  and  the  other 
in  accordance  with  Equation  3.36,  it  is  not  possible  to  state  which  of 
the  two  would  predict  more  accurately  what  the  heights  of  ripples  are 
in  laboratory  channels.  Again,  it  is  not  certain  that  a valid  formula 
can  be  derived  from  the  data  of  different  investigations  when  there  are 
differences  in  channel  geometry  and  in  granulometric  distribution  of 
grain  size.  There  are  two  reasons  for  the  difficulty.  One  is  that  the 
determinations  of  the  flow  quantities  and  the  dimensions  of  bed  features 
are  open  to  errors.  Secondly,  since  the  phenomenon  of  sand  transporta- 
tion is  a complex  phenomenon  the  parameters  entering  are  too  numerous 
and  the  error  of  measurements  prevents  the  evaluation  of  the  signif- 
icance of  a particular  single  parameter.  The  effects  of  granulometric 
distribution  are  not  well  known,  nor  are  the  effects  of  the  width  and 
form  of  channels.  In  the  views  of  Tanaka"5  the  width  of  channel  plays 
an  important  role,  and  its  significance  is  considered  in  a theoretical 
formulation  of  the  feature  geometry.  One  may  imagine  that  by  resorting 
to  theoretical  analysis  the  forms  of  the  proper  parameters  may  be  de- 
fined. But  this  may  not  always  be  possible.  For  example,  the  theory 
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of  Kennedy  predicts  that  L/H  should  decrease  with  (u/Uo)/U  . 

This  is  verified  in  the  Fort  Collins  data  obtained  with  sand  of 
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dj-g  = 0.93  mm.  In  the  same  data,  on  the  other  hand,  with  the  sand  of 
the  lesser  size  just  the  opposite  is  observed;  that  is,  L/Hq  increased 
with  (U  - UQ)/U  . Should  this  anomaly  be  explained  on  the  grounds  that 
d#  of  the  sands  lesser  than  d^Q  = 0.93  is  below  the  critical 
d#  = 16.5?  Even  when  tests  are  made  with  sands  of  dimensionless  size 
parameter  d#  less  than  l6,  we  sometimes  notice  differences  in  the 
behavior  of  the  size  parameters  as  reported  by  different  authors.  For 

4o 

example,  in  the  Shinohara  and  Tsubaki  results  with  sand  d = 0.21  mm 

2 

we  find  that  both  L/d  and  HQ/d  decrease  with  U /(s  - l)gd  , whereas 
the  opposite  was  true  in  the  Fort  Collins  results. 

77.  In  the  present  analysis,  it  was  not  possible  to  derive  the 
resistance  coefficients  of  alluvial  streams  in  the  form  established  by 
Alam  and  Kennedy  from  the  data  obtained  from  laboratory  channels.  In 
rivers  when  the  Froude  number  is  less  than  0.5,  A?  is  independent  of 
d/H  and  depends  on  U /gd  solely.  When  the  flow  Froude  number  is  less 
than  0.5  the  associated  bed  forms  consist  of  ripples  or  dunes.  In  con- 
trast with  the  above,  deductions  from  laboratory  testing  show  that  A? 
depends  also  on  d/H  in  addition  to  being  dependent  on  U /gd  . The 
analysis  involved  two  steps.  First,  it  was  sought  to  express  A in 
the  parameters  Hq/L  and  Hq/H  and  secondly,  to  relate  Hq/H  and 
Hq/L  separately  to  U^/gd.  It  may  be  of  interest  to  reexamine  the 
steps,  making  use  of  the  dimensionless  sand  size  d*  in  the  reductions. 
For  time  being  it  should  be  adequate  to  use  the  Alam-Kennedy  or  Einstein- 
Barbarossa  procedure  in  establishing  the  resistance  of  natural  rivers. 

The  Dune  Initiation  Criteria  in  Terms  of  Channel  Current 

78.  In  the  previous  development,  the  criteria  for  the  initiation 
of  ripples  and  dunes  were  expressed  in  terms  of  the  shear  velocity. 

For  the  ripples 

U#  = 0.72  [(s  - l)gv]1/3  (2.31  bis) 

and  for  the  dunes 
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with  the  restriction  that  d#  is  larger  than  3.2  but  smaller  than  16.5. 


In  the  Fort  Collins  data,  this  restriction  is  fulfilled  for  all  the 
tests  except  those  carried  out  with  sand  of  d^^  = 0.93  mm.  These  cri- 
teria can  be  readily  changed  to  new  relations  containing  U , the  cur- 
rent velocity,  instead  of  U#  . For  the  reduction,  it  is  sufficient  to 
obtain  the  ratio  U/U#  applicable  to  Fort  Collins  data.  An  explana- 
tory examination  showed  that  U/U#  would  vary  as  the  one-sixth  power  of 
H/k  . To  find  if  the  ratio  is  further  affected  by  the  Froude  number, 
the  graph  of  Figure  38  was  prepared.  The  quantity  U/U#  * (H/k)1^  is 


U/V  |H 

Figure  38.  Dependence  of  the  ratio  of  the  current  velocity  to  shear 
velocity  on  channel  depth  in  Fort  Collins  tests  (8- ft  flume) 

plotted  against  Froude  number  U/*^gH  . Data  are  differentiated  to  g^rin 
size  and  bed  features.  The  ordinate  of  each  point  represents  the  ave  \ge 
of  three  or  four  like  points  with  abscissa  values  very  close  to  eac’.; 
other.  It  is  obvious  that  the  graph  can  be  represented  as 


independent  of  bed  form  feature  and  Froude  number.  Substituting  in 
Equations  2.31  and  2.32,  the  criterion  for  the  initiation  of  the  ripples 
in  terms  of  current  velocity  is 

U = 2.59[(s  - l)gv]1/3(|)  ' (3.37) 

and  for  the  initiation  of  dunes 

U = 11.0  [ (s  - l)gd]1/5  (J)3/5  ' (3.38) 

Transitions  occur  at 

U = 3.1  C(s  - 1)]1/2  1/6  (3.39) 


a relation  that  follows  from  Equation  2.33.  By  dividing  the  two  sides 
of  Equation  3.38  by  VgH  and  putting  s = 2.67  , we  have 


U 

/9H 


(3. HO) 


This  allows  comparison  with  the  Vanoni  graphics  on  the  question  of  bed 
forms. The  above  relation  for  the  initiation  of  the  dunes  is  valid 
in  the  range  of  d*  greater  than  3.2  and  less  than  l6.5.  The  con- 
nection between  the  Vanoni  symbol  and  d#  is 

(s  - 1)R2  = d l (3. Hi) 

g 

Thus,  the  relation  in  Equation  3.H0  for  the  initiation  of  dunes  should  be 
valid  for  Rg  less  than  0.50  but  greater  than  H.H5.  Equation  3.H0 
requires  that  within  this  range  U/ VgH  should  vary  as  the  minus  one- 
third  power  of  H/d  . This  is  comparable  with  the  data  shown  in 
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Figures  2-5,  inclusive,  in  the  paper  by  Vanoni.  Although  in  these 
figures  the  paths  of  the  curves  separating  the  dunes  from  the  ripples 
are  somewhat  irregular,  on  the  average  U/t^gH  would  vary  as  the  minus 
one-third  power  of  H/d  . 

79*  Again,  in  the  similar  manner.  Equation  2.32  after  d is 
divided  by  v'gd  and  s is  put  equal  to  2.67,  for  ordinary  sands,  the 


critical  relation  for  the  initiation  of  dunes  would  be 


-0. 60 


(3.42) 


This  is  also  comparable  with  the  results  of  Vanoni"5  in  Figure  7 of  his 
paper.  The  equation  of  the  line  in  the  figure  separating  the  dunes  from 
the  ripples  is 


^==3.a(^d.0 


-0.75 


(3.43) 


The  agreement  between  these  two  results  is  close.  This  is  an  important 
matter,  since  these  equations  would  predict  what  the  effect  of  tempera- 
ture would  be  on  bed  forms.  It  is  to  be  remembered  that  the  validity  of 
these  relations  is  for  the  range  of  d#  greater  than  3.2  but  less  than 

16.5  or  for  the  range  of  R greater  than  4.45  but  less  than  50,  for 

6 

the  case  of  ordinary  sand,  s = 2.67  • For  another  comparison  of  the 

2 

results,  see  Figure  39  where  U#/(s  - l)gd  is  plotted  against  d*  . 

The  diamonds  (V)  are  reductions  from  the  straight  line  in  Figure  7 of 
Vanoni' s paper"5  separating  the  dunes  from  the  ripples;  the  full  circles 
are  reductions  from  the  Chabert  and  Chauvin  curve  in  Figure  17  of  the 
present  study  separating  the  dunes  from  the  ripples;  and  the  open  cir- 
cles are  the  reductions  from  the  Fort  Collins  data.  As  expected,  the 
points  from  the  Vanoni  study  align  themselves  on  a straight  line, 
whereas  the  alignment  of  the  points  from  the  Chabert  and  Chauvin  curve 
and  those  from  the  Fort  Collins  tests  is  curvilinear.  The  curve  that  is 
drawn  through  the  latter  points,  although  close  to  the  straight  line 
drawn,  suggests  that  the  effect  of  temperature  is  more  significant  for 
sands  with  d*  values  less  than  16.5  but  not  differing  much  from  it. 
Further  research  in  this  problem  would  be  desirable. 


PART  IV:  CONCLUSIONS 


80.  The  Manning's  ng  of  nuclear-created  channels  may  be  sup- 
posed to  consist  of  two  parts,  one  due  to  the  asperities  of  the  ejecta 
and  the  other  due  to  the  expansions  and  the  contractions  in  the  channel 
caused  by  the  disperse  explosions.  In  the  channels  of  basalt  ejecta, 
the  combined  ng  is  estimated  to  be  0.0762;  with  rhyolite,  0.0609  and 
with  clay  shale,  0.051*+  for  the  excavation  data  shown  in  Table  7. 

81.  In  a conventionally  created  channel  assuming  that  the  bed 
material  consists  of  Atlantic  muck  and  bed  features  are  absent,  the 
Manning's  n of  a wide  channel  60  ft  deep  would  be  0.013  for  a current 
velocity  of  10  ft/sec  and  0.015  for  a current  10  times  smaller.  In  the 
presence  of  bed  features,  the  resulting  combined  Manning's  n according 
to  the  method  of  Alam  and  Kennedy  would  be  0.015  for  a current  of 

10  ft/sec  and  0.030  for  a current  10  times  smaller.  . 

2 -*-/  -3 

82.  The  dimensionless  grain  parameter  d*  = [(s  - l)g/v  ] d 
introduced  by  Valembois  and  then  reconsidered  by  Bonnefille  plays  an 
important  role  in  the  sediment  motion.  Here,  s is  the  ratio  of  the 
density  of  the  grain  material  to  the  density  of  water;  g , the  constant 
of  gravity;  v , the  kinematic  viscosity  of  liquid;  and  d the  size  of 
the  grains,  the  median  diameter.  Fort  Collins  tests  indicate  that  when 
d#  is  greater  than  16.5,  the  initial  sand  entrainment  is  first  followed 
by  a flat  bed  and  then  by  dunes  as  the  current  velocity  is  increased 
without  the  intermediate  ripple  formation.  Where  d#  is  less  than  3.2, 
ripples  are  followed  by  transitions  without  the  intermediate  formation 
of  dunes.  In  the  range  3.2  < d#  < 16. 5 the  criterion  for  the  forma- 
tion of  dunes  is 


/(s  - 1 )gd 


= 3.06  d#°'9 


where 


U#  = the  critical  shear  velocity 
d = the  median  grain  diameter 

118 


i 


This  is  an  approximation  and  the  actual  fraction  showed  marked  deviations 
from  the  above.  Since  the  criterion  is  predictive  of  the  changes  pro- 
duced by  temperature,  these  deviations  need  to  be  studied  further. 

Putting  in  the  above  s = 2.6?  , the  resulting  expression  of  the  cri- 
terion agrees  reasonably  well  with  the  findings  of  Vanoni. 

83.  For  the  study  of  the  form  resistance  of  bed  features  in  terms 
of  their  characteristic  dimensions,  the  Sadron  formula  is  adopted.  In 
the  formula,  the  resistance  is  related  to  the  velocity  U prevailing 
at  the  apices  of  the  bed  forms,  ripple  or  dune.  Assuming  that  the  local 
velocities  in  terms  of  distance  from  the  bed  are  proportional  to  y^  , 
the  Sadron  formula  allows  for  the  form  resistance  coefficient  the 
expression 


where  A is  a constant  and  the  effect  of  B/H  is  neglected. 

84.  Fort  Collins  data  on  velocities  indicate  that  the  velocity 
index  m is  somewhat  higher  for  the  ripples  than  for  the  dunes.  There 
are  also  variations  with  the  grain  size.  Taking  averages,  m is  equal 
to  0.202  for  the  ripples  and  to  0.139  for  the  dunes.  For  the  studies 
made,  m is  taken  as  1.66,  an  average  from  all  the  determinations — 
ripples  and  dunes  alike.  As  regards  the  exponent  a , the  data  suggest 
that  it  is  permissible  to  take  a equal  to  unity.  Owing  to  scatter, 
this  selection  may  not  be  positive.  Adopting  a equal  to  unity,  factor 
A is  somewhat  larger  for  the  dunes  than  for  the  ripples. 

85.  When  the  crest  height  Hq  and  the  length  L of  the  bed 

features  are  expressed  in  terms  of  d , the  grain  median  diameter,  the 

2 

resulting  ratios  are  functions  of  (s  - l)U  /gd  unaffected  by  the  grain 

size  and  water  depth  ratio  d/H  . The  functional  dependencies  are  not 

the  same  for  ripples  and  for  dunes.  These  are  the  findings  from  the  Fort 

Collins  data.  There  is  an  indication  that  the  ratio  H /d  is  affected 

• o 

by  viscosity.  Owing  to  the  scatter  of  data,  the  resolution  is  not  def- 
inite and  for  the  analysis  made  this  possibility  is  ignored. 

86.  Introducing  the  depth  and  the  length  values  of  bed  features 
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in  terms  of  the  current  velocity  in  the  formula  of  resistance  involv- 
ing the  length  and  the  depth  of  the  features,  it  is  found  that  for 

laboratory  flumes  the  coefficient  depends  both  on  H/d  and 

2 ^ 

(s  - l)U  /gd  . In  rivers  as  found  by  Alam  and  Kennedy,  the  coefficient 

2 

is  a function  of  (s  - l)U  /gd  only  when  the  Froude  number  U/»/gH  is 
less  than  0.5*  As  the  Froude  number  of  the  laboratory  flume  considered 
was  also  under  0.5»  it  is  seen  that  the  resistance  processes  are  dif- 
ferent for  laboratory  flumes  and  natural  rivers. 

87.  Because  of  this  difference  in  the  behavior  of  laboratory 
channels  and  natural  rivers,  it  is  preferable  that  the  form  resistance 
of  deep  channels  carrying  fine  sediments  be  evaluated  by  the  Alam- 
Kennedy  method. 
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APPENDIX  A:  NOTATION 
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E’E1’E2’EL 


Domain  of  influence  of  the  frictional  action  of  the 
bottom,  A = A + A, 

OWl) 

Channel  cross  section 


Total  value  of  projected  area  of  all  the  asperities  over 
an  area  A 


Horizontal  projection  of  lee  faces  the  ripples 


Domain  of  influence  of  the  frictional  section  of  the  two 
walls 


Contact  area  between  current  and  isostatic  flow 


Area,  wetted  wall  surface  in  a crater  segment;  area  of 
separation  surface 


The  extreme  and  the  mid-cross  sections  in  a crater 
segment 


The  average  cross  section  of  the  i crater 


Channel  width 


Drag  coefficient;  velocity  at  a point  1 ft  from  bottom 
of  the  flume 


A coefficient  defining  A in  terms  of  H /L  and 
Hq/H  , Equation  3-3 


A coefficient  defining  A in  terms  of  L/H  , Equation 
3-10  ° 


A function  of  placement  of  the  shape  of  sands,  and 


Reynolds  number  U^d/v 


Median  diameter  of  sand  grains,  pebbles  or  rocks 


Bqnnefille  modification  of  Valembois  sand  parameter 

j _ n 


a;  = G 


Exposure  ratio 

Jet  dissipation  due  to  free  turbulence  in  a nuclear 
segment 

Darcy-Weisbach  resistance  coefficient,  f = ^A 


Darcy-Weisbach  resistance  coefficient  for  granular 
surfaces  of  bed  features,  f^  = UA1 

Darcy-Weisbach  resistance  coefficient  for  form 


resistance  of  bed  features,  fg  = UAg 


Resistive  force  of  a ripple  or  a dune  or  sand.  Also 
Froude  number  U/ ✓gH 


Resistive  force 
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Stress  Froude  number  F„  = 

Gravitational  constant 

2 3 

Valembois  grain  parameter  G = s - 1/v  gd  , s = p /p 

s 

Depth  of  water  in  channel;  also  range  of  tide 
Crest  height  of  ripples  or  dunes  measured  from  bottom 
Height  of  asperities;  Karman  constant 
Equivalent  sand  roughness 

A coefficient  in  the  resistance  equation  in  power  form 
for  granular  surfaces,  hydraulically  smooth 

Depths  of  channel  bottom  areas  of  uniform  asperities  k. 
and  of  Manning's  n^  ; i = 1,2, . . . , n 1 

Length  of  ripples  and  dunes;  usually  a length  of  channel 
segment 

Length  of  converging  and  diverging  parts  of  a crater 
segment  of  length  L . L = L^  + L^  • 

Constant  close  to  unity;  exponent  in  the  law  of  varia- 
tion of  separation  length  Jt.  ; exponent  in  the  power 
loss  of  velocity  distribution 

Manning's  roughness  coefficient 

Effective  Manning's  n for  a cratered  segment  due  to 
granular  surfaces  and  contractions  and  expansions  in 
the  segment 

Manning's  n for  composite  channel,  for  sidewalls  and 
for  bottom 

Wetted  perimenter  of  bottom,  p^  = B 

Wetted  perimeter  of  a complete  channel,  pQ  = B + 2H 

Wetted  perimeter  of  sidewall,  pw  = 2H 

Mass  rate  of  flow  of  sand  material  expressed  as  pounds 
per  hour  per  foot 

Volume  rate  of  flow  of  sand  material  expressed  as  cubic 
feet  per  hour  per  foot 

Discharge  through  a nuclear  segment 

Hydraulic  radius 

Hydraulic  radius  in  general 

Hydraulic  radius  corresponding  to  the  bottom  of  channel 
Hydraulic  radius  of  a composite  channel 
Hydraulic  radius  corresponding  to  the  sidewalls 


A2 


. r 


Grain  Reynolds'  number  U#d/v 

Hydraulic  radius  in  Einstein's  theory  associated  with 
the  granular  surfaces  in  a wide  channel 

Hydraulic  radius  in  Einstein's  theory  associated  with 
the  form  resistance  of  bed  features  in  a wide  channel, 

R = R1  + R2 

The  extreme  and  the  midhydraulic  radii  in  a crater 
segment 

The  average  hydraulic  radius  of  the  i^*1  crater 

Ratio  of  densities  of  grain  to  density  of  water, 
s = p /p  ; arc  length  measured  from  the  lowest  point  in 
the  section 

Energy  gradient;  energy  gradient  in  a composite  channel; 
energy  gradient  in  wide  channel 

Energy  gradient  in  Meyer-Peters*  theory  associated  with 
granular  surface  resistance  in  wide  channels 

Energy  gradient  in  Meyer-Peters'  theory  associated  with 
form  resistance  of  bed  features  in  a wide  channel 

Local  velocity  at  distance  y from  bottom 

Mean  current  velocity  in  channel 

Sand  entrainment  critical  velocity 

Current  velocity  at  the  apices  asperities 

Current  at  crests  of  bed  features 

Velocity  of  current  at  the  apices  of  sands 

Velocity  of  current  at  the  start  of  a nuclear  segment 

Shear  velocity,  U#  = A/p 

Shear  velocity  of  granular  surfaces,  = A-^/p 

Shear  velocity  of  bed  features,  U#2  = A2/p  , 

T = T1  + T2 

Fall  velocity  of  grains 
Horizontal  distance 

Vertical  distance  from  channel  bottom 
Depth  of  water  in  channel 

Average  length  of  separated  area  in  the  direction  of 
flow  in  a segment 

Length  of  separation  area  in  the  direction  of  flow 
observed  at  the  water  surface  in  an  expanding  segment 
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Angle  of  inclination 

Boundary  layer  thickness,  6 = 11.6V/U# 

Einstein's  symbol  for  effective  asperity.  Equations  2.38 
and  2.1+4 

Fall  of  energy  line 

Resistance  coefficient;  2A  = U2/U2  , 4l  = f 

Resistance  coefficients  of  composite  channel,  of  channel 
walls,  and  of  channel  bottom 

Resistance  coefficient  of  granular  surfaces, 

2XX  = (U#1/U)2 

Force  resistance  coefficient  of  bed  features, 

2A2  = (U*2/U)2 

Kinematic  viscosity  of  water 

Density  of  water 

Density  of  granular  material 

Gradation  number  2a  = d^g/d^g  + dg^/d^g 

Channel  stress 

Average  shear  of  bottom  surface,  = pgR^S 
Channel  average  stress  in  composite  channels 
Average  shear  of  sidewalls,  = pgR^S 
Stress  due  to  granular  surfaces,  = pgR^S^ 

Stress  due  to  bed  features,  t = pgRS  , 

S.  = S + S0  , S = S R b 

b 1 2 b o o 

Functional  symbols 

Average  wetted  parameter  in  a crater  expansion  segment 
Shields ' parameter 

Einstein's  transport  function  ip'  = 1.67d^,-/R^S 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Keulegan,  Garbis  Hovannes 

An  estimate  of  channel  roughness  of  interoceanic  canals  / 
by  Garbis  H.  Keulegan.  Vicksburg,  Miss.  : U.  S.  Waterways 
Experiment  Station  ; Springfield,  Va.  : available  from 
National  Technical  Information  Service,  1978. 

124,  4 p.  : ill.  ; 27  cm.  (Technical  report  - U.  S-  Army 
Engineer  Waterways  Experiment  Station  ; H-78-13) 

Prepared  for  U.  S.  Army  Engineer  District,  Jacksonville, 
Jacksonville,  Florida. 

References:  p.  121-124. 
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